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HAUPTAUFSÄTZE 
Über rechteckige Tragflügel sehr großer Tiefe. 


Von P. Kriesis in Athen. 


Die Luftkräfte auf einen unendlichtiefen Tragflügel können aus der Wider- 
standszahl der zweidimensionalen ebenen Platte berechnet werden. Bei 


Seitenverhältnissen ergibt sich dann ein einfacher Ausdruck 


für den Auftrieb, der gut mit entsprechenden Messungen von Winter 
übereinstimmt. 


1. Der rechteckige Tragflügel mit kleinem Seitenverhältnis ist in der letzten Zeit von 
W. Bollay') und K. Wieghardt?) behandelt worden. In beiden Untersuchungen werden 
die Luftkräfte aus den gebundenen Tragflügelwirbeln, welche senkrecht auf der Anströmungs- 
geschwindigkeit V stehen, und dem Kutta-Joukowskischen Satz berechnet. Die ab- 
gehenden Wirbel liegen bei Wieghardt in der Richtung von V, während Bollay annimmt, 
daß sie mit dem Tragflügel einen bestimmbaren Winkel bilden. Wieghardt hat mit den- 
selben Ansätzen den Fall untersucht, daß die Hinterkante im Unendlichen liegt, während die 
Saugkante im Endlichen bleibt (a. a.O. S. 266), und Bollay hat den Grenzfall eines zu beiden 
Seiten unendlichtiefen Tragflügels behandelt (a.a.O., S. 33 bis 35). 

Wir werden im folgenden zeigen, daß das Problem des unendlichtiefen Tragflügels nicht 
mit den Methoden der Tragflügeltheorie behandelt werden darf, da es unmittelbar mit der 
zweidimensionalen Strömung um eine ‘ebene Platte zusammenhängt. Es ergibt sich dann für 
die Normalkraft auf den beiderseitig unendlichtiefen Tragflügel ein Wert, welcher mit den 
experimentellen Ergebnissen übereinstimmt, und aus diesem 
ein Näherungswert für den Auftrieb eines Tragflügels end- 
licher, aber sehr großer Tiefe bei kleinem Anstellwinkel. 

2. Wir betrachten in Bild 1 einen zu beiden Seiten un- 
endlichtiefen Tragflügel, in dessen Ebene die Achse Oy senk- 
recht auf der Anströmungsgeschwindigkeit V steht, Oz parallel 
zu den Seitenkanten in der Mitte der Spannweite b, und Ox 
senkrecht auf Oy und Oz ist. Wenn dieser Tragflügel unter 
dem Anstellwinkel a mit der Geschwindigkeit V angeströmt 
wird, so muß man für die Berechnung der Luftkräfte die 
Druckverteilung und das Geschwindigkeitsfeld als Funktionen 
des Ortes und der Zeit bestimmen. Die strenge Lösung dieses 
Problems genügt, als instationäre Strömung aufgefaßt, den 
Navier-Stokesschen Differentialgleichungen: 


) W. Bollay: 7. angew. Math. Mech. Bd. 19 (1939), S. 21. 
°, K. Wieghardt: Z. angew. Math, Mech. Bd. 19 (1939), S. 257. 
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der Kontinuitätsgleichung: + + 0, und den Randbedingungen: =Pp,, =Vsina, 
te —=0, m»—Vcosa im Unendlichen, und: u=vr=mw=0 auf dem Tragflügel, also auf der 


Fläche: - o<2<+tm.,. 


Unser Problem wird aber wegen: 
in folgende zwei getrennte Probleme aufgespaltet: 


1. Es sollen p (x, 9, u (x, 9,8), v(x,9,#) aus den Navier-Stokesschen Differential- 
gleichungen in zwei Dimensionen und der entsprechenden Kontinuitätsgleichung, die wir nicht 
wieder hinschreiben, berechnet werden, mit den Randbedingungen: p& =p,, v» =0, u» —=Vsina 


im Unendlichen, und v=v=0 auf der Linie: — sys 


Die Lösung dieses Problems ist aber gerade‘ die zweidimensionale Strömung um eine 
ebene Platte der Breite b, die senkrecht mit der Geschwindigkeit Vsin«a angeströmt wird. 


2. Es soll die Funktion w (x, y, t£) berechnet werden, welche der Differentialgleichung: 


genügt, und die Randbedingungen: w=0 auf der Linie: -3 sys und =Vcosa 


erfüllt. Dabei müssen die Funktionen u (x, y,t) und v (x,9,t) aus dem Problem 1 eingesetzt 
werden. 


3. Man erkennt auf diese Weise, daß erhebliche Geschwindigkeitskomponenten in der 
2-Richtung auftreten, welche sogar örtlich und zeitlich rasch veränderlich sein werden und 
bei kleinen Anstellwinkeln viel größer als die anderen Komponenten sind. Dieser Umstand 
macht jedoch für die Druckverteilung absolut nichts aus. Es genügt also für die Berechnung 
der Druckkräfte auf den Tragflügel nur die Druckverteilung auf das Profil einer ebenen Platte 
bei senkrechter Anströmung zu kennen, welche ja experimentell gut bekannt ist. 

Die Reibungskräfte auf der Tragflügeloberfläche sind dagegen von den Reibungskräften 
auf die ebene Platte im zweidimensionalen Problem verschieden, da jetzt in dem Tensor der 
” und 7 auftreten, die erst berechnet 
werden können, wenn das Problem 2 gelöst ist. Bei allen diesen Strömungen kann man aber 
bekanntlich die Reibungskräfte auf die Körperoberfläche gegen die Druckkräfte vernachlässigen. 
Auf der Tragflügeloberfläche wirken daher dieselben Luftkräfte wie auf einer ebenen Platte, 
deren Profil in der &-y-Ebene senkrecht mit der Geschwindigkeit u» =Vsina angeströmt 
wird. Die gesamte Druckkraft ist also eine Normalkraft, die in der x-Richtung senkrecht auf 


Deformationsgeschwindigkeiten auch die Ableitungen 


die Tragfläche steht und pro Längeneinheit der Flügelsehne den Betrag hat: N zug b. 


Hierbei ist o die Luftdichte und c,„, die Widerstandszahl, die bekanntlich für die ebene Platte 
bei senkrechter Anströmung den Wert 2 hat’). Wegen uo=Vsina folgt hieraus für den 
Koeffizienten der Normalkraft: 


V’b 

Diesen Wert haben aber gerade die Experimente von Winter*) nach Flügel’) für einen 

unendlichtiefen Tragflügel ergeben. Unsere Ableitung läßt erkennen, daß er für alle Anstell- 

3) Vgl. z. B. O. Flachsbart: Z. angew. Math. Mech. Bd. 15 (1935), S. 32. 


4, H. Winter: Forsch. Ing.-Wes. Bd. 6 (1935), B, S. 40. 
5) G. Flügel: Schiffban, Bd. 30 (1929), S. 336. 
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winkel von O bis z gilt. Man sieht ferner, daß der Kutta-Joukowskische Satz in dieser 


Strömung nicht gilt, da sie abgerissen ist, und daß die Luftkräfte durch das Auftreten der 
Wirbel in der z-Richtung hervorgerufen werden, während die Wirbel in der y-Richtung eine 
ganz nebensächliche Rolle spielen. Die Annahmen von Bollay sind also in der Strömung 
um den beiderseitig unendlichtiefen Tragflügel nicht erfüllt. Seine Rechnung führt übrigens 
in diesem Fall, trotzdem sie den richtigen Wert cy=2sin?’a ergibt, auf einen Widerspruch 
betreffend die Gestalt der abgehenden Wirbellinien, auf welchen Bollay auch hinweist. 


4. Weniger einfach als der vorhin behandelte Grenzfall ist der Fall eines einseitig un- 
endlichtiefen Tragflügels.. Man muß dabei unterscheiden, ob die Vorderkante oder ob die 
Hinterkante im Endlichen liegt. 

Den ersten Fall hat Wieghardt für kleine Anstellwinkel mit den Methoden der Trag- 
flügeltheorie behandelt. Seine Untersuchung ergibt für den Auftrieb des gesamten unendlich- 
tiefen Tragflügels den endlichen Wert: A=0,789 0 V’b?’sina. Dieses nicht sehr wahrschein- 
liche Ergebnis rührt davon her, daß nur an der Vorderkante gebundene Wirbel in der 
y-Richtung auftreten, während weiter hinten die Normalgeschwindigkeit — Vsina auf der 
Tragflügeloberfläche allein von den abgehenden Wirbeln induziert wird. Nach dem Kutta- 
Joukowskischen Satz sollten daher etwas hinter der Vorderkante keine Luftkräfte mehr 
auftreten. Das trifft aber in Wirklichkeit nicht zu; denn die starken abgehenden Wirbel der 
Vorderkante lösen sich sehr bald auf, und es tritt an ihrer Stelle allmählich die abgerissene 


Strömung um den beiderseitig unendlichtiefen Tragflügel auf, die bei 2z=+» voll aus- 


gebildet ist. 

Zwischen dem Gebiet, wo diese Strömung herrscht und der Vorderkante wird also nicht 
mehr die Normalkraft N=sin’«aoV?b pro Längeneinheit der Flügelsehne wirken, sondern 
es werden erheblich stärkere Kräfte auftreten. Wir fassen alle diese zusammen, und nennen 
AT den Auftrieb, der sich insgesamt an der Vorderkante der Normalkraft N pro Längeneinheit 
der Flügelsehne überlagert. Dieser zusätzliche Auftrieb rührt von den gebundenen Wirbeln 
in der y-Richtung her, die nach Wieghardt den Beitrag AY,= 0,7890 V?’b’sin«a liefern, 
und von den hohen Unterdrücken, welche die starken, abgehenden Wirbel der Vorderkante 
bei ihrer Auflösung an den Seitenrändern des Übergangsgebietes hervorrufen. Wir nehmen 
an, daß diese zusätzlichen Kräfte der abgehenden Wirbel des Übergangsgebietes linear von 
der Intensität der gebundenen Wirbel abhängen, und daher auch dem Auftrieb A; propor- 
tional sind. Dann ist: 

wo £ eine Konstante bedeutet, die wir später bestimmen werden. 

Da der Tragflügel unendlich tief ist, und auf ihn außer dem Auftrieb A? auch die 
Normalkraft N=sin’aoV?b pro Längeneinheit der Flügelsehne wirkt, erfährt der Tragflügel 
also tatsächlich einen unendlichen Gesamtauftrieb. 


5. In dem Fall, wo die Hinterkante im Endlichen liegt, werden in ihrer Umgebung die 
gebundenen Wirbel in der y-Richtung nur bei großen Anstellwinkeln eine Rolle spielen. Bei 
kleinen Anstellwinkeln müssen wir aber annehmen, daß zu Beginn der Bewegung ein Anfahr- 
wirbel abgespült wird, so daß von den beiden Seitenkanten zwei Wirbelzöpfe in Strömungs- 
richtung abgehen und in der Nähe der Hinterkante nur ganz schwache zusätzliche Wirbel in 
der y-Richtung auftreten. Es wird also in diesem Falle keine wesentliche Änderung der 
Druckverteilung entstehen, und man wird die Luftkräfte wie beim beiderseitig unendlichtiefen 
Tragflügel berechnen können. Der Unterschied mit dem im vorigen Abschnitt -behandelten 
Fall liegt gerade darin, daß dort jenseits des Tragflügels keine Wirbel in der z-Richtung 
vorhanden sind, so daß an der Vorderkante starke Wirbel in der y- Richtung auftreten müssen. 

6. Die Strömung um einen Tragflügel endlicher, aber sehr großer Tiefe wird bei kleinen 
Anstellwinkeln von den vorhin behandelten nicht sehr verschieden sein. Wir nehmen wieder 
an, daß sich der Normalkraft N=sin’aoYV?b pro Längeneinheit der Flügelsehne ein Auf- 
trieb A, überlagert, welcher nur von den gebundenen Wirbeln in der y-Richtung abhängt, die 
sich an der Vorderkante ausbilden, so daß für alle Seitenverhältnisse, für welche A, = 4% 
ist, auch Ar= APR gilt. Den Anteil A, von Apr kann man aber für das jeweilige Seiten- 
verhältnis direkt aus Bild 8 der Untersuchung von Wieghardt entnehmen (a. a. O., S. 268). 


Man findet, daß für alle Seitenverhältnisse 4 die kleiner als 5 sind, Ag gleich dem asymp- 
totischen Grenzwert = 0,789 0 V? sin a gesetzt werden kann. Für diese Seitenverhältnisse 
ist also die gesamte Normalkraft auf den Tragflügel wegen (1): 

N=Af cosa+sin’aoV?’bt 


sinacosa+tsin’aoV’bt. 
1* 
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Und der Koeffizient der Normalkraft: 


=2 sin? a+ 0,7897 sin 2a. 
—zoV?bt 
2 
wo die Konstante £ noch bestimmt werden muß. & 
Wenn man cy über 3 aufträgt, so stellt bei > 
festem { die Gl. (2) eine Schar von Geraden dar, „,, 
deren jede einem gewissen Anstellwinkei ent- 
spricht. In Bild 2 sind diese Geradenscharen für „’ 
gestrichelt und für {=2,15 ausgezogen ge- 
zeichnet worden. Man sieht, daß die von Winter I —_ 
im Bereich 0 < 7 < gemessenen Werte von cx, 
besonders bei den kleineren Anstellwinkeln, mit ot 
guter Näherung auf der ausgezogenen Geraden | 
schar liegen, so daß £=2,15 gesetzt werden muß. > 
Auf diese Weise wird also die Beziehung 7. | 7 
=0) —=2sin’a und die von Gl. (2) geforderte | __—7__+-” 
lineare Abhängigkeit des Koeffizienten cy von L------ 
# 
sin 2a und Im Intervall durch die 
Versuchsergebnisse von Winter gut bestätigt. * 
Für größere Seitenverhältnisse darf aber Gl. (2) 
nicht mehr angewandt werden, denn es ist dann Die Geraden der GI, (2) mit t=3.15. 
Ag AG — - - - Die Geraden der G]. (2) mit {=1. 


sina“ wie Bild 8 der Unter- 
suchung von Wieghardt erkennen läßt. Tatsächlich fallen auch die gemessenen Werte 


1 
von cy für ?> 3 unter die ausgezogenen Geraden. 


Für Seitenverhältnisse 0 < 7< - können wir jetzt auch den Auftrieb des Tragflügels 


berechnen. Es ist wegen {=2,15: 
A= +cosasin’aoV?bt 
Und für kleine Anstellwinkel «: 


Hieraus folgt die Auftriebszahl und der Auftriebsbeiwert: 


43407, 
_ b 


Bei den Tragflügeln sehr großer Tiefe hängt also ze linear von dem Anstellwinkel « und 


dem Seitenverhältnis 3 ab*®). 524 


Eingegangen am 17. 11. 1943. 


6), Die in Anmerkung ®) und 5) angeführten Arbeiten von Winter und Flügel sind dem Verfasser in Athen 
nieht zugänglich gewesen. — Die in Bild 2 angegebenen Versuchsergebnisse von Winter sind aus der Arbeit von 
Bollay (a.a.0O., S.31, Fig. 11) entnommen worden. 
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Biegungsschwingungen eines in beiden Enden unterstützten, 
punktförmig belasteten Balkens. 
Von T. G. Owe Berg in Bofors (Schweden). 


Es werden die Schwingungen eines Balkens, die durch Gl. (1) dargestellt 
werden können, unter Einwirkung einer plötzlich angelegten punktförmigen 
Masse M, die mit der Geschwindigkeit v,—=0 bzw. v,7£0 antrifft und an 
dem Balken haftet, untersucht. Die zwei niedrigsten Schwingungszahlen 


werden berechnet und als Funktionen von % und S graphisch dargestellt. 


Wir betrachten im folgenden einen in beiden Enden 2<=0 und x=L unterstützten 
Balken, der von einer im Punkte = plötzlich angelegten Masse M zum Schwingen ange- 
regt wird. Die Differentialgleichung der Schwingung ist, wenn wir von der statischen Durch- 
biegung unter Einwirkung des Eigengewichts des Balkens absehen, 


mw 
mit den Randbedingungen 
=0 für ISO. . .(2a), . . . ab. 
Da nur für ist, gilt für alle die einfache Schwingungsgleichung 


Für £=£ erhalten wir nach Multiplikation von (1) mit dx und Integration zwischen den 
Grenzen E—e und 


E+e 
3 3 
Wir lassen jetzt e>0 und setzen lim mas Dann wird für = 
> 


Wir teilen den Balken in zwei Gebiete ein: L.O<r<!; Il. Ezx-=L, und führen 
zwei Funktionen w, und w,, eine für jedes Gebiet, ein. Im Punkte «= gelten die Stetig- 
keitsbedingungen 

und die Bedingung 


wer 


Eine partikulare Lösung von (3) ist 
w„—=[4A„Coshayi (n=1,2), 


wobei der Kürze wegen a=/ 84 gesetzt wurde. Die Randbedingungen (2b) und (6) geben 


Sinhayi(L— 
sinayi(L- 


®sinayALSinh ay 


B, = 


D,=— 


| 
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CoshayiL 
:SinhayiL' 
sinayAs 
cosayAL sinayÄE 
’sinayALSinhayis 


Die allgemeine Lösung kann daher geschrieben werden 


D,=A 


+® 
sinay4LSinhayi(L— &Sinh ayix — 


+» 

| sinay A 


. (Sb), 


wobei die Funktion F(4) aus (5a) bestimmt werden soll. Da der Nenner der Integranden 
in Gl. (8a) und (8b) nur reelle Nullstellen hat, kann die Integration längs der reellen 
Achse durch Integration längs des in Bild 1 gezeigten geschlossenen Weges I, ersetzt werden. 
Der Weg L setzt sich aus einer Parallelen zur reellen Achse, unmittelbar unter dieser, und 
aus einem Halbkreis, dessen Halbmesser R ist, zusammen. Lassen wir R>x, so gehen die 
Integranden von (8a) und (8b) gegen Null. 
Es wird also 
|=|=2ni Res., 
wobei Z’Res. die Summe der Residuen bezeichnet. 
Wegen 


erhalten wir aus (5a) 


sinay A LSinhay A(L — &)Sinh ayAs— 


Sinhayd: E EI sinay L-Sinh ay AL:Sinhayas 
_Mg \ 


7 dd. 
Es wird also ä 
1 Mg 
Fü) = | 
sin ay AL Siohayi Yilsin a yALSinha y 


-SinhayA&—SinhayiLsina yA(L— &)sina y 


Wir führen die Transformation 
a L V = 


ein und verändern den Integrationsweg nach Bild 2. Dann wird 


Bild 1 (links). Bild 2 (rechts). 


B,= w; 
wo 
wo 
(83), 
: 
und 
Das 
i= 
aufl 
1 für 
an 
0 
1 
R 
6 
. (10). 10 
(10) 
0 
1 
| L’ 2 
6 
10 
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i 
n,()=5, 


Isin pSinh p Sinh pZ — Sinhpsinp 


p* IsinpSinhp Sinhp Sinhp — Sinhp sinp ginp 


wo M,=oFL die Masse des Balkens ist. Für w,(x,f) gilt ein entsprechender Ausdruck, 
wo im Zähler x und £ vertauscht werden sollen. 


Die Integration ergibt (nur das halbe Residuum im Pole p=0 soll mitgerechnet werden) 


EI d M co IL? 
dp „sin Sinhp„+ 5M, Pn |sin p„ Sinh —— L Sinh 


— Sinh p„ sin Sinpn 


und entsprechendes für w, («,f). Dabei sind p,„ die reellen Nullstellen des Nenners 
Sinhp + 57; sin p Sinh p Sinhp  — Sinhpsinp" |. (13) 
sınp Sin sinpSinh p Sinhp 7 inhpsinp Sinp7, (13) 


In Gl. (12) ist das erste Glied der rechten Seite der Beitrag des Poles p=0 zum Integral. 
Das zweite Glied stellt die Summe der Beiträge der Pole p„ #0 dar. Setzen wir in Gl. (12) 
{=0 ein, so erhalten wir w,(x,0)=0, da sich die beiden Glieder der rechten Seite dann 
aufheben. Da die Reihe gliedweise nach f deriviert werden darf, ist auch offenbar 0 


für =0. Die Lösung (12) genügt also den Anfangsbedingungen. 


Zahlentafel 1. Zahlentafel 2. 
Pa 


0,05 | 0,10 | 0,15 | 0,20 | ı 0,25 | 0,40 
0,95 0,0 | | 05 | 00 | 0 | 000 


| | 
| 20 | | 20 | | 
5,981 15,474 — 5,191 
5,708 5,084 4,914 4,957 5,105 | 
5.055 4587| — | — 


4,447 14,343 | | 


4.066 , 4,927 


a 


DD 


a8 


Pi 


| | 9 | 
0,50 0,05 | 0,10 | 0,15 020 | | 0,25 | 0,40 | 0,50 
| 


0,75 | | , 0.90 | 0,90 | 0,85 | 0,80 | 0,75 | 0,60 | 


| | | 
8,994 | 6,849 5,683 : 35,77 30,51) — 12695 27,91 135,24 
8,289 5,523 | 4,393 32,58 25,17 124,15 24 5726,06 — | 4m 
6,381 | 3,531 | 2,719 — — | - 
5,331 | 2,816 | 2,140 261982 — — 


7 
an, 
), 
(12) 
), 
n 
d 
e 
0,10 | 0,25 ' 0,08 | 
| 0,90 | 0,70 | 950 | | 0,97 | j 
| MIMo | 
n a|ı 0 an | 2n | 
1 | a | 2,999 | 2,617 | 2,384 1221 — | 
2 | | 2,879 | 2,350 | 2,096 2 — | 
6 | 2,526 | 1,879 | 1,649 6 — | 
| 2,309 | 1,678 | 1,463 10 | 2272| — | — | 
20 | 2,005 1423 | 1.937 | — | - | | 
oo 4.730 | 4,010 — 1486| — 
pi 
0,10 | 
| 0 | 
| 
0 n?| 
1 
2 | m? | 
6 | 
10 
20 
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Bild 3. 


u Zahlentafel 1 gibt die be- 
rechneten Werte von p, und p»} 


für verschiedene Werte von 
6 M, 


und LT: Die entsprechenden Kur- 


für p? sind in Bild 3 und 4 
| ven für pt sind in Bi un 
dargestellt. Zahlentafel 2 und 
3 Bild 5 und 6 geben Entsprechen- 
des für p,, p} bzw. p}. 
} es Pr pi zw Pi | 
Die Reihe (12) konvergiert 


sehr schnell. Mit einem Fehler 
kleiner als 1°) können wir statt 
(12) schreiben 


Bild 4. 


Das Resultat (12) bzw. (14) 


kann leicht für eine mit der Ge- 
fi | schwindigkeit v, auffallende Masse 
2 erweitert werden. Wir nehmen 
/ dabei an, der Stoß sei vollkommen 
\\ unelastisch, und setzen statt (2) 


dt 


für z=E£, t=0. 


| | Um dieselben Anfangsbedingungen 
wie früher zu erhalten, führen wir 
die Funktionen 


ein. Diesen Funktionen sollen die Bedingungen (2) und (6) genügen. Dagegen gilt statt (5a) 


(5b). 
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Bande Berg, Biegungsschwingungen eines punktförmig belasteten Balkens 


Die Funktionen «, und uw, 
wurden so gewählt, daß EN-050 — 


du, 
9 


für t=0 und 


0 


Sie geben Funktionen w, und w,, 
deren Ableitungen 


dm, 
7 


für i=0, z=E. 


Diese Funktionen w erfüllen nicht 
etwa eine Bedingung 


Sie geben eine bestimmte Ge- 
schwindigkeitsverteilung zur Zeit 
t=0, und zwar eine mit der 
statischenDurchbiegungsverteilung 
identische Geschwindigkeitsver- 
teilung‘). Physikalisch bedeutet 
diese Anfangsgeschwindigkeitsver- 
teilung eine unendlich große 
Ausbreitungsgeschwindigkeit der 
Störung. Dieselbe Eigentümlich- 
keit zeigt die Lösung (12), was 
aus Gl. (14) besonders deutlich 
hervorgeht?). Entsprechend (9) 
setzen wir 


0 


für 


Im Zähler von (10) fügen wir demgemäß . i 1 en Fu hinzu und in (12) 

sin Sinh p, Sinh p, Sinh p„sinp,„ in PT 
sin 27, 


Pa dp. Isin PnSinhp„+ \sin p„Sinhp,„ Sinh 


— Sinh p„sin p„ sin 
bzw. in (14) 


aI®L 
s 


pi a? 2 


Eingegangen am 15. 9. 1943. 
1) Handbuch d. Phys. VI, Berlin 1928, S. 550. 
2) Flügge: Z. angew. Math. Mech. Bd. 22 (1942), S. 313. 
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Der Spannungszustand der biegungssteifen Kegelschale mit 
linear veränderlicher Wandstärke unter beliebiger Belastung. 
Von H. Nollau in Dresden. 


Die vorliegende Arbei‘ entwickelt eine Methode zur Berechnung der biegungs- 
steifen Kegelschale bei beliebiger Belastung und enthält die Berechnung der 
Kreiszylinderschale und der Kreisplatte als Sonderfälle. Die Darstellungs- 
weise und die Erläuterung der Theorie durch ein ausführliches Zahlen- 
beispiel bezwecken, die Ergebnisse der Abhandlung technisch verwendbar 
zu machen. 


I. Einleitung. 


Die Literatur enthält bisher vollständige Angaben über den Spannungs- und Verschiebungs- 
zustand der Kegelschale nur für drehsymmetrische Belastung. Unsymmetrische Lastfälle 
wurden nach der Membrantheorie behandelt. 

Die Aufgabe dieser Arbeit besteht in der Berechnung der biegungssteifen Kegelschale 
für beliebig auftretende Belastung der Ränder bei kräftefreier Schalenoberfläche. 

Die Integration der Differentialgleichungen des Problems gelingt jedoch nur bei Kegel- 
schalen, deren Wandstärke zum Abstande von der Kegelspitze proportional ist. Diese Ein- 
schränkung hat für Kegelschalen mit gleichbleibender Wandstärke keine Bedeutung, da der 
Membranspannungszustand durch vorgeschriebene Randbedingungen nur in einem schmalen 
Bereiche verändert wird. Dasselbe trifft für Rotationsschalen zu, deren Meridian mit ge- 
nügender Genauigkeit durch eine Gerade ersetzt werden kann. 

Vergleichende Berechnungen sind mit der Kegelschale konstanter Wandstärke #, durch- 
geführt worden. Es hat sich gezeigt, daß sich bei den verschiedenen Randbedingungen Ersatz- 
schalen mit {=t’-y derselben Mittelfläche angeben lassen, deren Extremwerte M,„ als 
wesentlicher Schnittkraft mit denen der Kegelschale konstanter Wandstärke übereinstimmen, 
Die Betrachtungen galten nur für rotationssymmetrische Belastung, können aber auch auf 
höhere Harmonische übertragen werden. Die Ersatzschale hat bei sehr 
dünner Wandung nahezu dieselbe Randstärke wie die Schale mit t,, da- 
gegen schneiden sich die Leibungen bei dicken Schalen nach Bild 1 in 
einem Abstande m vom Auflager. Dabei ist die Strecke m bei Belastung 

unteren oben 
des Randes nach 
richtet sich nach den Randbedingungen. 

a) Gelenkige, unverschiebliche Auflagerung. 

Bei Vernachlässigung der Querdehnung 4 können als Randbedingungen 
der biegesteifen Schale verwendet werden: 

= — Ngyembran M, 
In diesem Falle ist etwa 
bis 
Dabei bedeutet Z die Länge der Randwelle der Abklingungsfunktion für 
die nullte Harmonische der Kegelschale mit ’=t,/l. Wird die Quer- 
dehnung berücksichtigt, so sind folgende Randbedingungen zu erfüllen: 


abzutragen. Die Länge der Strecke m 


Ep — Ey Membran» M, 
Sollen die Ringdehnungen der beiden Schalen am unteren Rande gleich 
sein, so ist die Randstärke der Ersatzschale geringer anzunehmen, als dies 
bei Übereinstimmung der Randringkräfte der Fall sein müßte, da der im 
Vergleich zu f, vergrößerten Randstärke eine größere Randringkraft ent- 
spricht als — Npyempran- Die umgekehrten Überlegungen gelten für den 
oberen Rand. Der Abstand m kann jetzt etwa mit 


L 
angenommen werden. 
b) Starre Einspannung. 


€&9 = — Membran» — MW Membran} 
oder bei u. —=0: 


‚ 
Ny Ny Membran» = — W YMembran- 
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Die Wahl der Ersatzschale richtet sich, abgesehen von der Wandstärke f,, nach dem 
Verhältnis e,/w’ oder N,/w’ und damit nach einer Funktion der äußeren Kräfte. Allgemeine 
Angaben sind also ausgeschlossen. Indessen liefert die Vergleichsschale, deren Wandstärke 

unteren kleineren 

oberen größeren 
als max M, der Schale gleichbleibender Wandstärke ,. Das Ergebnis befriedigt, wenn als 
Stärke der Ersatzschale #, im Abstande 


Rande mit t, übereinstimmt, für M, einen nur unwesentlich Grenzwert 


m=( bis 


gewählt wird. 

Die Schnittkräfte M,, N,, Q, infolge der Randzwängung besitzen gegenüber M,, Ny 
keine Bedeutung für den Spannungszustand. 

Ist an beiden Rändern mit Störungen zu rechnen, so sind zwei verschiedene Ersatz- 
schalen für die Abklingung vom unteren Rande und für die Abklingung vom oberen Rande 
aufzusuchen. 

Die Rechnung wird unter den folgenden Voraussetzungen durchgeführt: 

Daher können die Normalspannungen o, und die Gleitungen yyz=Yyy; gleich Null gesetzt 
werden. Die Punkte einer ausgezeichneten Schalennormalen liegen infolgedessen nach der 
Formänderung wiederum auf einer Normalen zur verformten Mittelfläche. 

2. uvu,n<t. 
Diese Annahme bringt es mit sich, daß man bei den Verschiebungsgrößen nur lineare Glieder 
zu berücksichtigen braucht. 


Bild 3. Die Spannungen am Schalenelement. 


Bild 2 (links). 


Zur Lösung der Aufgabe werden die folgenden 
Bezeichnungen gebraucht (Bild 2): 


a halber Öffnungswinkel des Kegels. 
r Hauptkrümmungsradius der Schalenmittelfläche. 
a Abstand eines Punktes J des Schalenmittelfläche von der Drehachse. 
Schalenstärke. 
!' Winkel zwischen der Innen- und Außenleibung der Kegelschale. 
9,9=»cosa,z Koordinaten der Schale. 
y Abstand des Punktes J von der Kegelspitze. 
® Breitenwinkel des Punktes J. 
p Projektion des Winkels 9 auf den Mantel des Normalkegels zur Schalenmittelfläche. 
Positiver Drehsinn von p und » entgegen dem Uhrzeiger. 
z Koordinate in Richtung der Schalennormalen mit der Mittelfläche als Ursprung und 
dem positiven Richtungssinn entgegen zum Krümmungsmittelpunkt. 
u,v,w Verschiebungen des Punktes J der Schalenmittelfläche in Richtung der Koordinaten y, 9, 2. 
Bild 3 zeigt die positiven Richtungen der Spannungen, Bild 4 die positiven Richtungen 
der Schnittkräfte. 
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Die Differentialopera- 
tionen werden folgender- 
maßen abgekürzt: 


f 
of of 
usw 


Bild 4. 


II. Aufstellung der Differentialgleichungen. 


Das Spannungsproblem der Schalen umfaßt die Berechnung von 10 Schnittkräften und 
3 Verschiebungen, sodaß 13 Bedingungsgleichungen zur eindeutigen Berechnung der Unbe- 
kannten nötig sind. Das Gleichgewicht der inneren und äußeren Kräfte am differentialen 
Abschnitt liefert mit 6 Bedingungen die 6 Gleichungen (2a—f). Hierzu treten die 8 Gleichungen 
des Verschiebungsansatzes der Schnittkräfte (8a—h), in den die geometrische Verträglichkeit 
und das Hookesche Gesetz eingehen. Es können also zusammen 14 Beziehungen zur Be- 
rechnung der 13 Unbekannten aufgestellt werden. Der Überschuß von 1 Gleichung ist jedoch 
nur scheinbar, denn die Gleichgewichtsbedingung (2f) erweist sich als Identität, wenn man 
die Schnittkräfte durch die Verschiebungen ausdrückt. Dies kommt von der paarweisen 
Gleichheit der Schubspannungen 7,9 =7,,, die ja aus einer ganz entsprechenden Momenten- 
gleichung gefolgert wird. Es bleiben also 13 unabhängige Beziehungen, die sich zu 3 homo- 
genen Differentialgleichungen (9a--c) mit den Unbekannten u,», w als Grundlage der Theorie 
zusammenfassen lassen. 


A. Gleichgewichtsbedingungen. 


Die Spannungen werden durch Integration über die Schalenstärke zu Schnittkräften zu- 
sammengefaßt, die auf die Länge 1 der Mittelfläche in den Schnitten y= const und = const 
bezogen werden: 


+t/2 
Längskräfte: N,= N= 
12 
+t/2 +t/2 
Schubkräfte: = | \wyda; 
+t/2 +t12 
Querkräfte: Q,=— \r,,(1+zir)dz, %=— \1.da; . (la—k). 
—t/2 
Biegemomente: M=-— | oyzdz; 
-i1 — 12 
+t12 +t12 
Drillungsmomente: =— \ Myy=— | wy2dz | 
—i12 


Die Koordinate z wird also nicht gegen den Krümmungshälbmesser r vernachlässigt, 
und damit nicht von vornherein lineare Spannungsverteilung angenommen. Der allgemeine 
Ansatz bereitet keine Schwierigkeiten, so daß die Gleichgewichtsbedingungen vollständig be- 
friedigt werden können. 
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Nach Bild 4 lassen sich die folgenden 6 Bedingungen für das Gleichgewicht der Schnitt- 
kräfte am differentialen Abschnitt der Kegelschale anschreiben: 


1. Gleichgewicht der Kräfte in Richtung der Meridiantangente 
yNy+Ny+ N;yeota— N, =0 
2. Gleichgewicht der Kräfte in Richtung der Breitenkreistangente 
N; Nyp+Nyo+ Npy— 
3. Gleichgewicht der Kräfte in Richtung der Schalennormalen 
eota+ &cota+yQ,+ Q,=V 
4. Gleichgewicht der Momente um die Meridiantangente 
M,cota+yMyga + + My y—-y = 0 
5. Gleichgewicht der Momente um die Breitenkreistangente 
Mr —yQ,=0 
6. Gleichgewicht der Momente um die Schalennormale 
My „eota+y(Nyp — )=0 


B. Elastizitätsgesetz und geometrische Verträglichkeit. 


Nach den Annahmen in Abschnitt I bleiben von den 6 Komponenten des Formänderungs- 
zustandes des elastischen Körpers nur drei: Die Meridiandehnung e,, die Ringdehnung e, und 
die Gleitung yyp=yy, übrig. Sie sind durch das Elastizitätsgesetz mit den Spannungen o,, 
0 und verknüpft: 


E 
&y), 


E 


Typ Yuy 


Um die Schnittkräfte als Funktionen der Verschiebungen u, v, w der Mittelfläche an- 
zuschreiben, wird die Verzerrung des Flächendifferentials durch den Verschiebungszustand 
der Mittelfläche ausgedrückt. Dies geschieht in 2 Schritten, von denen der erste die Ver- 
schiebungen der Mittelfläche mit der Verzerrung des Differentials der Mittelfläche verbindet, 
der zweite mit der Verzerrung beliebiger Flächendifferentiale der Schale. 


1. Verzerrungsgrößen der Mittelfläche. 


B 

(4a—e). 
u ‚vw 

| 


2. Verzerrungsgrößen &,, &y, Yyy im Querschnitt. 


Die Verzerrungsgrößen einer Mantelfläche vom Krümmungsradius r+z gehen aus den 
Gl. (4a-c) hervor. Dabei wird r durch r+z und «, v, w durch «,, v,, w, ersetzt: 
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Breitenkreiselement 


D 
BZ 


Bild 5. Bild 6. 


Der Unterschied zwischen 
den Verschiebungen «;, v., 
eines Schalenpunktes A im 
Querschnitt und den Verschie- 
bungen des zugeordneten Punk- 
tes A, in der Mittelfläche ent- 
steht nach Bild 7 aus der Be- 


wegung der Schalennormalen 
Y AA, bei der Verformung 
r+z w 
Merichanschnät (6 a-c). 
W; — 
Bild 7. (m + 


Die Beziehungen zwischen e,, &, y,y und u, v, w lassen sich nunmehr unter Benutzung 
der Gl. (5) und (6) folgendermaßen anschreiben: 


[23 
r+z 
r 


z 
(utana + ztana 2, 


3. Schnittkräfte. 
Die Verbindung der Gl. (7) mit dem Hookeschen Gesetz (3) liefert die Spannungen 


und mit den Ansätzen (la—d, g—k) die Schnittkräfte. Die Integrationen nach z werden in 
den logarithmischen Gliedern nach t/r entwickelt und nach der 3. Potenz abgebrochen. 


als Dehnungssteifigkeit und 
u) 


als Biegungssteifigkeit der Schale lassen sich 8 Schnittkräfte als Funktionen der Verschie- 
bungen u, v, w in folgender Form ausdrücken: 


+ 


eota+ cota) eota, 
Yy Yy 


N, = y cota+ y cot a +K cot a+ al, 


. 1 .r . (Sa-d), 
y eota+v K eot’a a+ a eot?’ a 


1l—u 


Ney="5"D 


w „[w" 2 w 2 w 3 
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M=K|n 


cota— 
Yy 


+ cot’a+ eotat um“ 
(&e—h). 


„[w'’ w 


„|w’ u cot?a , vcota v’ cota 
Mu K |" ota— 


Dieser Ansatz gelingt jedoch nicht für Q, und @,, da yy; und yy, nach den grund- 
sätzlichen Voraussetzungen Null gesetzt worden sind. 


C. Differentialgleichungen des Spannungsproblems der Kegelschale. 


Der Ansatz (8) enthält keine Angaben über die Querkräfte, die daher aus den 5 Gleich- 
gewichtsbedingungen (2a—e) für die Schnittkräfte entfernt werden. Dazu eliminiert man sie 
in den Gl. (2d) und (2e) und setzt sie in die Gl. (2b) und (2c) ein. Gl. (2a) kann in der ur- 
sprünglichen Form angeschrieben werden. 


YNy+ Ny+ Nyyeota— N, 


N,cota+yNyy+ Nyy — — My, cota — = -cota cot a=0, 


 (2a—c)'. 
Nyeota— Mi 2 y My +" cota+ 


Diese Gleichungen liefern mit (8) und mit k=K/D=1t?/12 drei Bedingungen für das 

Gleichgewicht der Schnittkräfte als homogene simultane Differentialgleichungen 4. Ordnung 
mit «, », w als veränderlichen Größen. 


+ 
3—u 3 
+y? "cot’a — 5 “cota— u — wcota kw cot’a— km’ cota (9a), 
3 2 3 
cot’a— cot’a — kucot’a— kweot’a—0 


kw cot’a+ 


1— 1 
+yp 


1 


(9b), 
+Yy — u) kw cot’a (1 — u) kv’ cot? a) 
3 3 
+ + uw’ cota]+ 
kv” cot?a— kw’ + cot’a+ucota+ mcot’a 
3 
Ru km” yrcotta— ku’ cota+ko (9e). 


+km"cotta+ ku" cot’a+2kw” (cot*a + 2cot? a) — 
3+u 


cot?a+ ku(2cota + cot’a) + km (cot!a +2 cot’a) = 0 


Die Abhängigkeit der Wandstärke von y geht in den Faktor k ein. Drei Fälle sind 
von Bedeutung: 


\, 
w; 
im 
ık- 
en 
| 
. | 
| 


Z. angew. Math. Mech. 


16 Nollau, Spannungszustand der biegungssteifen Kegelschale Band 24 Nr. 1/1944 
1. t=const : k= 15 const; 
2 


Die mathematische Lösung für die Ansätze 1 und 2 mit konstanter und beliebig linear 
veränderlicher Wandstärke bereitet Schwierigkeiten, da sich die Integration mit bekannten 
Funktionen nicht vollziehen läßt und Reihenansätze sehr langwierige Rechnungen erfordern 
würden. Dagegen gelingt für die Kegelschale, deren Wandstärke nach 3. proportional dem 
Abstand von der Spitze ist, eine geschlossene Lösung. 

Ansatz (8) verwandelt sich mit t=!’y, 


in 
N,=Df[yw + cota + cota)] — K y’ w” cota, 
N=D[u+v 


cota+yv' Klymw” cot’a— cot’a+vcot’a— yv’cot’a], 


a— w' cot’a+ w'cot’a], (l1a-h), 
M,= — y? u’ cota — uyv'cot? a], 
K[ymw" cot?a + ywceot’a+ y’w’ + yucota+uy'w”], 


und das System (9a—c) erhält nach Division durch %? die folgende Form: 
cot a] — 


+ 


2 


v 'cot’a+(u—k)w’ cot a] — 


u cot’a(l+kcot’a)+(l — Boot cota+ 


En“ u + w cota—0 


y eu +3K& cot?a) > 


+v" cot’a— (1 — cot’a)e +(1+3(1— u) k) w' cot’«—=0 


ky w’’+ y k[-u”cota+8 + 


+ cot a+2 mw” cot’a+ (11 +34) + 


+(2— u) ku” cot’a+(1+k (cot’a — 1)) ueota+(1 (3 — 5u)) v’ cot’a+ 


+ kw“ kw” cot’a+(1+ k (cot? — 1)) weot’a=0 

Die Arbeit von A. Pflüger [1]') über die Stabilität der Kegelschale enthält ähnliche 
Gedankengänge, die dort ausführlich dargestellt worden sind. Der Differentialansatz (13) 
seiner Untersuchung läßt sich nach Streichung der Lastglieder in die Gl. (12) überführen. 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. 
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III. Auflösung der Differentialgleichungen für die Wandstärke t= dy. 


A. Charakteristische Gleichung. 

In jedem Term der Gl. (12) ist der Grad der Ableitung nach % gleich dem Grade der 
Potenz der Veränderlichen y. Es sind also Eulersche Differentialgleichungen, die einen 
Potenzansatz in y ermöglichen, der die Auflösung des Systems in Verbindung mit einem 
Fourieransatz in @ zuläßt. 


n 
> Un sin. p | (13a-e), = By (14a—c), 
| 
> ? 
J 
n=0,1,2,3,...; 


Die Fourierschen Reihen (13) werden mit n/cosa angesetzt, da der Winkel p bei 
einer Umdrehung des Meridians von 0 auf 2x nur um 2rcosa anwächst. Die Substitution 
i=#'+1/2 empfiehlt sich für die einfachste Berechnung von 4. 

Die Differentialgleichungen (12) ergeben mit den Ansätzen (13) und (14) für jede Har- 
monische n die 3 linearen Gl. (15a—c), in denen die Konstanten A}, B;, C} und der Ex- 
ponent A, als Unbekannte auftreten. Sie werden als Matrix angeschrieben. Da sich die 
folgenden Rechnungen immer auf eine bestimmte Harmonische n beziehen, kann der Index n 


wieder wegfallen. 
B* 


(15a—c). 


In 
|, 


4 2 sin’a sin?’a 


—=kcota 


In 7—5 
|, 
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Das System ist homogen und A*, B*, C* sind deshalb nur dann von Null verschieden, wenn 
die aus seinen 9 Koeffizienten d,x gebildete Determinante verschwindet. Sie führt auf eine 


symmetrische Gleichung 8. Grades für A +4. Die in k quadratischen Glieder sind klein 
und werden vernachlässigt. Auf diese Weise läßt sich folgende Bestimmungsgleichung für A 


anschreiben. 


Die Beiwerte g, bis g, sind nach Multiplikation mit er ! 


»=-[7—-4u 
+ 


+ 84) cwt’a— (I — 8 u) cott a +2 (2 al’ + 
+[13 —124+2(9— 11 u) cot?’a — (3+8 u) cot a — 2(4 — u) + 


+4(1 + cot? a)? (16 a-e). 


9,8 


- - +24 42) — + 124°) a+ 


a+ jgeotta+ 


u)eot* a 
+ 1212) a+[ 3.) a— alm'+ 


—4u+(7— 12 u) cot?a+3(1 — 4 u) cot*a — (3 +4 u) cot*a — 2 cot’ 


+ (1+ cot?a)* | 


Die Gl. (16) enthält nur gerade Potenzen in A und ist daher eine Gleichung 4. Grades 
für 4°. Die Wurzeln [2] können reell und komplex sein und in folgender Gliederung auftreten: 


1. alle 8 Wurzeln sind reell; 
2. alle 8 Wurzeln sind komplex; 
3. 4 Wurzeln sind reell, 4 komplex. 


So entsteht der folgende allgemeine Ansatz für die Exponenten A,, in dem bei einem reellen 
Wurzelpaar Ä,,,,, der N,=0 und zu Setzen ist: 
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B. Allgemeine Lösung des Ansatzes. 
1. Abklingungsfunktion. 
a) Mathematische Lösung. 
Die allgemeine Lösung der totalen Differentialgleichungen jeder Harmonischen n lautet: 


r—=8 


r=8 
r=1i 
r—=8 
3 
r=]1 


Da die Exponenten 4, reell oder konjugiert komplex sind, gilt dasselbe für die Konstanten 4A*, 
B*, C*. Damit jedoch in den Ausdrücken für die Verschiebungen «„, ®,, ®„ nur reelle 
Glieder vorkommen, wird der Ansatz (18) in den komplexen Gliedern r, r+ 1 umgeformt. 


(r=1,3,5,7), 


1 1 
= 2 + Cr, 
r 


> 
= y  ?[C%(eos nr 


+ +, (cos (n, In y) —i sin In y))], 


+ C7 4) cos (nr In y) + — sin (nr In y)] 


0, — C,,) . . . . (19). 
C, und C,;, sind reell, also C} und C},, konjugiert komplex. 


[C,cos(nriny)+ Cr; ,sinyriny) (r=1,3,5,7%).. . . (20a). 


sı 


1 
 ?[C,cos(n, In y)+ €, sin (n, In y)]+ 


1 

eos (n, In y)+ C, sin (niny)]+ 

+Yy 2[C, cos (n, In y) + C, sin (n, In y)] + 


2 [C,cos (n,Iny) + C, sin (n, In y)]. 


Dieselben Gleichungen gelten für «„ und v„ mit den Integrationskonstanten A, und B,. 
Wenn Ä, reell ist, können an Stelle von 20 (a—c) unmittelbar die Gl. (18) benutzt werden. 


2. Beziehungen zwischen den Integrationskonstanten. 


Zu jeder der 8 Lösungen 4,, r=1-—-8, der Gl. (16) gehören nach (14) die drei Inte- 
grationskonstanten A}, B}, C}. Sie sind durch die Matrix (15) derart miteinander verknüpft, 
daß eine gewählt werden kann und damit die beiden anderen bestimmt sind. Die acht 
Integrationskonstanten C}, (r—=1— 8), der Verschiebung w als wichtigster Komponente werden 
gewählt und daher A}, B} mit den ersten beiden Gl. (15) als Funktion von C/ angegeben. 
Die dritte Gleichung kann zur Nachprüfung verwendet werden. Die Beiwerte der Gl. (15a, b) 
sind einfacher als die von (l15c), so daß die numerische Auswertung auf diese Weise weniger 
Fehlerquellen enthält. 


In jeder Lösung r wird: 
2° 


2 
| 
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b,=dscdga— 

=d,4dpp — dpa 

Bei zunzenn Abklingungszahlen ist außerdem noch das Verhältnis der se Konstanten 

A,, B,, ©, der Lösung (20) notwendig. Die Beiwerte von (15a) und (15b) für 4,,,4,=&+in- 

sind: 


2 


2 n? 
1 ‚1+ 


sına 


Das obere Vorzeichen gilt für ET das untere für 4,+,=&,—in,. Die komplexen 
Koeffizienten d;x ergeben die Verhältniszahlen a,,r+,, Br,r+1, Er,r+. als letzte komplexe 
Größen der Rechnung. 
Denn nun erfolgt mit 


und (21) der Übergang auf die gesuchten reellen Beziehungen: 
Artı = Cr tar Cr 
Darin bedeuten die reellen Faktoren a und ß: 
(24). 


Arlr+ı br 1 —dr+1Cr 
2 Cr Cr4; 2 
Die numerische Auswertung wird erleichtert, wenn die konjugiert komplexen Zahlen a, 
b, c zur Bestimmung von a und f in ihren Real- und Imaginärteil aufgespalten werden: 
=aRtayi, 
b,,r+,=br + by 
Damit gehen die Gl. (24) über in: 


Ar +1 = 


arcr+taycy brer+byc, 
a D Pr= ? 
ch+c 


(25). 
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Die Abklingungsfunktionen «4, und v„ lassen sich dann für komplexe Lösungen 4, , + , folgender- 
maßen anschreiben: 


r—@r+, Cr4,) cos (nr Iny) +(ar+, Cr+arCr+,) sin (n,Iny)], 


(20b,c 
[dr Pr +, Cr+,) cos (n-Iny)+(Pr+, Cr + Pr + ,) sin (n-Iny)] 


C. Komponenten des Spannungs- und Verschiebungszustandes. 


Der vollständige Ansatz (11) für die Schnittkräfte enthält zahlreiche Glieder, die für 
das Ergebnis keine Bedeutung besitzen und in Zukunft vernachlässigt werden. Damit ent- 
stehen nur Fehler in der Größenordnung von 0,1 bis 1 v.H. 

Der verkürzte Ansatz (11) lautet g 


weota)] — Ky’mw”cot’a, 


M,=Kl[y’w" + u(ym” cot’a+y’ w')], 
+uy’w”], 
My, = Myy = (1— u) K[y’ w’ cot a— ymw' cot a] 
Die Qwerkräfte folgen aus (2d,e) mit 
—Mp+Miycota), 
(26, h). 


Hierzu treten bei einem kräftefreien EM die beiden Kirchhoffschen Ersatzkräfte nach 
Abschnitt III D: 


cot a, 
(%6i, k). 
T,=Q, cot a 


Die Schnittkräfte werden nach (26), die Verzerrungen der Mittelfläche nach (4) aus (20) ab- 
geleitet. Von den Verschiebungen wird außer «, v, m meist nur noch w’ gebraucht. 

Für die Zahlenrechnung werden dimensionslose Koordinaten o=yıl, mit l„= y,„ oder 
o=y/l, mit ,=y, verwendet, je nachdem Randstörungen untersucht werden, die vom 
unteren oder oberen Rande einer offenen Kegelschale ausgehen. 


Aus diesem Grunde werden die Ansätze (20) und (26) mit 
y=ol, 
auf die dimensionslose Koordinate o umgestellt. Die beim Differentiieren entstehenden Aus- 
drücke für die Komponenten des Spannungs- und Verschiebungszustandes werden als 
Funktionen von o und # in Zahlentafel 1 vereinigt. 
Die Beziehungen zwischen den Integrationskonstanten ©, für benannte und (, für un- 
benannte Koordinaten sind: 


C,=1-#(C, In) — C,+,sin 
sin (n, In )) 
(C, cos (n, In D)+C,+,sin(n, In 
(Cr +, —C, sin(n-Ind) 


(27); 
oder 
(28). 


D. Randbedingungen der genauen und der genäherten Lösung. 


Die allgemeine Lösung der randbelasteten Kegelschale enthält 8 Integrationskonstanten 
zur Erfüllung von 4 Bedingungen an jedem Rande. Sie werden also aus einem System von 
8 linearen Gleichungen berechnet. Der Spannungs- und Verschiebungszustand der geschlossenen 
Kegelschale ist durch die Integrationskonstanten C/ — C/ bestimmt. Die Integrationskonstanten 
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C;—C% sind Null, damit die Verschiebungen und Schnittkräfte an der Kegelspitze ver- 
schwinden. Ausnahmen sind bei n=0 und n=1 möglich (vgl. Abschnitt IV B). 

Die Abklingungszahlen A,,, und A,,, sind nach Abschnitt IVB im Regelfalle, verglichen 
mit A,,, und 4,,,, klein. Sie erzeugen Längsspannungen und dehnungslose Verbiegungen, die 
für n>2 nicht mehr in reiner Form auftreten. Die Konstanten der Anteile 3,4 und 7,8 
werden ebenfalls verhältnismäßig klein, so daß der Spannungs- und Verschiebungszustand in 
Randnähe mit großer Genauigkeit durch die großen Lösungsanteile 1,2 und 5,6 beschrieben 
wird. Der Abschnitt V4 enthält den rechnerischen Nachweis. Der Unterschied in der Dämpfung 
der beiden einem Rande zugeordneten Anteile kehrt die Größenverhältnisse jedoch um, so daß 
die kleinen Anteile der Lösung im mittleren Schalenbereich überwiegen können und dort die 
Membranspannungen nur unerheblich verändern. Aus diesen Gründen genügen häufig die 
Anteile 1,2, so daß für jeden Rand nur 2 Konstanten bestimmt werden müssen. Das sind 
C},, und C;,,, mit denen dann die wesentlichen Randbedingungen erfüllt werden. 

Der Verzicht auf die schwach gedämpften Anteile der Lösung führt zu einer weiteren 
Vereinfachung der Rechnung. Das Verhältnis (l, — l,)/t„ ist bei den im Bauwesen verwendeten 
Schalen oft so groß, daß sich/obere und untere Randstörung nicht mehr beeinflussen. In 
diesem Falle werden C},, und C;,, unabhängig voneinander berechnet. Die Brauchbarkeit 
einer Näherungslösung läßt sich abschätzen. Hierzu dient außer (l, — l,)/t, die Randwellen- 
länge L der Abklingungsfunktion der höchsten vorkommenden Harmonischen. Die Rand- 
störungen gelten als unabhängig, wenn (l,—1,)/t„ >50 oder wenn (,—1)>L-2L. Die 
Amplituden sind dann auf weniger als 1/500 ihrer Randwerte abgeklungen. 

Die Erfüllung der statischen Randbedingungen bei einem kräftefreien Rande gelingt 
nicht ohne weiteres, da 4 Integrationskonstanten in die 5 Bedingungsgleichungen 

N,=N, ... . (29) 
eingehen würden. Diese Schwierigkeit ist durch G. Kirchhoff beseitigt. dust, der die 
3 Schnittkräfte N,y, Q, und M,y nach (26i,k) zur Ersatzschubkraft S,, und Ersatzquer- 
kraft T, zusammenfaßte. Die Kräftegruppe S,y, T, ist zur Kräftegruppe N,y, Qy, Myy äqui- 
valent und hat nur in einer Zone von der Breite £ eine andere Wirkung als diese. 

Daher gilt nunmehr für den kräftefreien Rand: 


N,=S8Sy =Ty=M,=0 . (30). 


IV. Diskussion des Ergebnisses. 
A. Angenäherte Berechnung der Wurzeln der charakteristischen Gl. (16). 

Die Einführung von Zahlen in die charakteristische Gl. (16) zeigt, daß oft ihre strenge 
Auflösung entbehrt werden kann. 

Die Querzahl u besitzt erfahrungsgemäß keinen großen Einfluß. Es genügt daher, sie 
in den wichtigsten Anteilen der Koeffizienten g, bis g, zu berücksichtigen, nämlich in den 
mit der großen Zahl 1/k multiplizierten Summanden. 

Die Beiwerte g,, 9,, 9. Sind außerdem viel größer als g,, so daß die Wurzeln A an- 
genähert aus den folgenden zwei quadratischen Gleichungen für Ä? berechnet werden können: 


Is 
32 I _ (31 a, b). 
Die Gl. (31a) liefert die vier großen Wurzeln A,,,, A,,., die Gl. (31b) die vier kleinen 
Wurzeln A,,,, A,,s- Sie erfüllen die charakteristische Gl. (16) um so genauer, je mehr das 
Produkt der Gl. (31) mit dieser übereinstimmt. 
Die Diskriminante der Gl. (31a) für die großen Eigenwerte 


ist, abgesehen von Ausnahmefällen bei Kegeln mit sehr großem Öffnungswinkel, negativ und 
liefert dann komplexe Abklingungszahlen. Dagegen wird die Diskriminante der Gl. (31b) für 
die kleinen Eigenwerte 
(33) 
je nach der Belastungsfunktion (n) und den geometrischen Verhältnissen (k, a) kleiner oder 
größer als Null und löst Störungsfunktionen mit oder ohne Vorzeichenwechsel aus. 
Die Kurven in Bild 8 beschreiben bei den Harmonischen n=2 bis n=5 die Abhängig- 
keit zwischen a und f’ für die Bedingung 4,—=0 mit «=0 und damit den Übergang von der 
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aperiodischen zur periodischen Schwingung der vier kleinen Glieder 3, 4 und 7,8 der Ab- 
klingungsfunktion. Die graphischen Beziehungen zwischen a und f’ bei vorgeschriebenem n 
unterteilen also die Kegelschalen mit reellen und komplexen Wurzeln 4,,,,;,s. Kegelschalen 
mit komplexen Wurzeln sind Punkte im oberen Bereiche zugeordnet 

Die Angaben im Abschnitt IV B mit den Bildern 9 bis 12 zeigen den Grad der An- 
näherung für A mit (31a,b) an die strenge Lösung für A mit (16). Die Näherung ist danach 
bei den niedrigen Harmonischen, mittleren Öffnungs- und kleinen Leibungswinkeln gut. 
Sie versagt bei flachen Kegelschalen und wird nach Bild 12 um so besser, je mehr der 
Kegelstumpf in den Zylinder übergeht. Das Kriterium 4,=0 für die reelle oder komplexe 
Form der kleinen Wurzeln A,,,,-,, ist also bei den plattenähnlichen Kegelschalen unbrauch- 
bar, während es mit wachsendem Basiswinkel immer schärfer wird. 


| ——— -Niherungslösung 


& 


5 
N 
T 
+ 4 - A 
cor? & Bild 9. Abhängigkeit der Abklingungszahlen 
Bild 8. Die Form der kleinen Lösungen auf Grund der Be- von den Glieuern n der harmonischen Reihe. 
dingung Ss a—=450%, 


B. Abhängigkeit zwischen der Abklingung der Randstörungen und den 
Konstanten n, a, f', «a der Lösung. 

Der Einfluß der Ordnung n der Belastung und der geometrischen und elastischen Eigen- 
schaften a, t', u der Kegelschale auf die Randstörungen wird durch Schaulinien geklärt. Ihre 
charakteristischen Eigenschaften bleiben, wie sich durch zahlreiche Untersuchungen gezeigt 
hat, bei anderen Konstanten der Kegelschale erhalten. 


1. Einfluß der Lastwechselzahl n. 

Bild 9 zeigt die Funktion Ä=f(n). Die Sonderfälle n=0 und n=1 führen zu den 
reellen Wurzeln A,,,=+1/2 und 4,,,=+3/2 der charakteristischen Gl. (16) und damit zu 
den folgenden Exponenten 4;,,, 4,,, in (18): 


+12 | -ı2 | 


Die Exponenten =0 und 4 kennzeichnen verformuugslose Bewegungen [3]. 4 —=—1 
und 4=—2 sind die Exponenten des Ansatzes (14) für singuläre, in der Kegelspitze an- 


greifende Einzelkräfte P und Punktmomente M. Die Konstanten 07, €, sind dann allein von 
der Intensität dieser Kräfte abhängig und können gegebenenfalls aus den Gleichgewichts- 
bedingungen bestimmt werden. 

a) Rotationssymmetrie des Spannungs- und Verschiebungszustandes mit n =. 

Der Exponent A/—=0 bedeutet eine Verschiebung der Schale in Richtung der Kegelachse mit 


v=0, w=-— A,tana, 
der Exponent 4, = + eine Drehung um die Kegelachse mit 
. v=By, w=®. 
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i,;=—1 beschreibt die Abklingung bei einer Einzelkraft P in der Drehachse, , = —2 die 
Abklingung bei einem Punktmoment M um die Drehachse. Sie erzeugen Spannungszustände 
ohne Biegung. Die Lösungen mit 4, und 4, lassen sich nicht durch den Fourieransatz (13) 
angeben, sondern werden erst nach Vertauschung der Funktionen sin und cos erhalten. 
Division der vollständigen Hauptgleichung (16) durch die vier. reellen Eigenwerte liefert 

die Gl. (34) zur Bestimmung der Abklingungszahlen A,,, und 4,,, für die Biegung am oberen 
und unteren Rande: 

j 9 2 2 1- uw 2 65 2 7 2 

Die Wurzeln sind, falls nicht ganz flache Kegelschalen untersucht werden, komplex mit 
und Dabei ist 


und . (85). 


9 
oe=T—u(2+cota), cot?’a +7 —Tu+3u? 


Diese Lösung ist zuerst von E. Meißner [4] auf anderem Wege entwickelt worden. 
Sie erscheint dort gekürzt wegen der Loveschen Vereinfachungen [5], auf die in dieser Arbeit 
verzichtet worden ist. Die Ergebnisse unterscheiden sich nur in Gliedern von untergeordneter 
Bedeutung. 

Bei der Berechnung der Randstörungen kann der Lösungsanteil 3 mit C},4, meist ver- 
nachlässigt werden. Er ist als formänderungslose Bewegung ohne Einfluß auf den Spannungs- 
zustand und bei der großen Steifigkeit der Schale gegen Längskräfte von geringem Einfluß 
auf den Verschiebungszustand. Der Anteil 4 scheidet aus, die Anteile 7, 8 können als Mem- 
branzustände unabhängig behandelt werden. Die Konstanten C},, und C?, sind daher Null. 
Zur mathematischen Formulierung der Randbedingungen dienen also lediglich C},, für den 
unteren und (5,, für den oberen Rand, die mit den Abklingungszahlen 4,,,,,,.=+(E+tin) 
die Störungsfunktion bestimmen. 

Das Verhältnis A}/C;=a,/c, wird mit der Vereinfachung 

(r=1,2,5,6) 
berechnet. 

b)n=1. 

Der Exponent 4, =0 bedeutet eine Verschiebung winkelrecht zur Kegelschale, parallel 
zur Meridianebene mit 

u=—Bsinaa v=B}, mw=-—Bjcosa, 
während A,=+ 1 die Drehung um die Normale durch die Kegelspitze zur Meridianebene 9=0 
beschreibt. Hierbei ist 

u=0, v=—(cosa-y, 
4,=—1l und ,=—2 führen bei einer waagerechten Einzelkraft P in der Nullmeridianebene 
und einem Kräftepaar M mit waagerechtem Vektor in der W°-Meridianebene zu Membran- 
spannungen. Die übrigen 4 Eigenwerte sind in folgender Gleichung enthalten: 
17 273 


+9- gut 19u+3 = 0 (86). 


Real- und Imaginärteil ihrer Lösungen 4,,,=+£f+in, ),,,=—+in sind: 


17 


(37). 


1— u? 17 273 
Die Anteile 3, 4, 7, 8 der Lösung scheiden ebenso wie bei n=0 für die Randbedingung aus, 
so daß deren Berechnung wiederum allein mit den Anteilen 1,2 und 5, 6 erfolgt. 


2. Einfluß, des Leibungswinkels !'=)12%k und der Querzahl u. 
Bild 10 zeigt die Funktion Ä=f(f). Der Einfluß der Querzalıl « läßt sich mit großer 
Genauigkeit nach Abschnitt IV A durch eine Abänderung von k beschreiben. w besitzt durch 


2 
die Verkoppelung SE mit & denselben Einfluß auf die Abklingungszahlen wie f'. 
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3. Einfluß des Kegelöffnungswinkels 2a nach Bild 11. 
Alle Eigenwerte erhalten bei flachen Kegelschalen eine Größenordnung, so daß die 
Näherungsansätze (31) versagen. 


C. Grenzübergänge gegen Platte und Zylinder. 


Kreisplatte mit =’: y und Zylinderschale mit t=const sind Sonderfälle des Problems 
und daher in dem allgemeinen Ansatz enthalten. Die Eigenwerte für die Kreisplatte ent- 
stehen mit cota=0 bei beliebigem n aus (16), für n=0 aus (35) und für n=1 aus (37). 


Der Übergang von der Platte über die Kegelschale zum Zylinder im Sinne von Bild 12 
schließt den Potenzansatz (14) aus, da a gegen Null geht, so daß ein Grenzübergang auf den 
Exponentialansatz e* der Lösung für den Zylinder nötig ist. Er lautet mit Bild 13: 


t 
a \cote]ik tan« 


Mit 


cot« 


entsteht der Zylinderansatz: 


lim =w,= (,- eXla-iz, 
Daher ist 


Der Grenzübergang auf den Zylinder bildet die Kontrolle der Hauptgleichung (16) und 
der Sonderlösungen (35) und (37). Die Gl. (16) für die Kegelschale geht mit 


lim A, =4,cota und lim k,=k,tan’a, k,=#?/12a?, 
a>0 


in die folgende, ihr zugeordnete N für die Zylinderschale über: 


— 2(2n’— —D)+1 40) 


—2n’[2n — (4 — u)n’+(2— u)]i2+n*(n? —1)’=0 
liefert nach (35) die Eigenwerte +&,+in, mit 


n=1 liefert nach (37) 


Die Ergebnisse stimmen bis auf die sehr kleinen zu IHRE Glieder mit den 
Lösungen der Spannungstheorie der Zylinderschale [6] überein. 


D. Wellenlänge und Dämpfung der periodischen Abklingungsfunktion. 
Der Verlauf einer beliebigen Schnittgröße S über den Meridian nach Zahlentafel 1 ist 


Ss=J o’r (fr ©, — fr + Cr + ) cos (n, In o) + (fr+ + fr + ,) sin (nrIn o)] (8). 


r=1,3,5,7 


Die Funktion läßt sich mit 6G,siny,, @,cosy, umformen: 


[fr sin (nr Ino+ yr) —fr+ , (n,Ine+ yr)] 


Die Summe der beiden trigonometrischen Funktionen läßt sich nun durch die Phasen- 
verschiebung x, in eine reine Sinusfunktion verwandeln: 


(44). 
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Bild 10. Abhängigkeit der Abklingungszahlen 4 vom Wandungswinkel ft’ und von der 


Querzahl u(n=31,a = 45°). 
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Bild 11. Abhängigkeit der Abklingungszahlen A vom Kegel- 
öffnungswinkel 2a, 


n=3, u—0. 
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Bild 12. Verlauf der Abklingungszahlen beim Übergang vom 
Zylinder zur Platte. 
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dabei ist 
__fr+: 


Hieraus entstehen für die Wellenlänge L, und das Verhältnis v, zweier aufeinanderfolgender 
Amplituden einer einzelnen Schwingung r einfache Ausdrücke. 


Bezeichnen 94 > ei die Orte zweier aufeinanderfolgender Nullpunkte. oder Ampli- 
9, m» Co,m+tı 


tuden eines Anteils r der Abklingungsfunktion vom rande, und damit 
die halben Wellenlängen, so ist bei einer u 
Abklingung vom Unterrande (r=1,3): 


Ou,m 
Ou ‚m+i 
Abklingung vom Oberrande (r=5,7): 


Qu,m — 


So.mtı _ lo,m+1— 00, m —1). 
00,m 


Die Wellenlänge ist demnach proportional zum Abstande von der Kegelspitze: 


Lu=21,(Ou,m 


Sie beträgt bei einer 
Abklingung vom Unterrande (r=1,3): = 2 ym (1 — ""r), 
Abklingung vom Öberrande (r=5,7): L,=21,0,m (""—1) =2 ym (er —1) 

Das Verhältnis », zweier aufeinanderfolgender Amplituden (verschiedenen Vorzeichens) ist 

„Pr 
„=e " (47) 
mit pr nach Zahlentafel 1. 


V. Spannungszustand in einem kegelförmigen Kühlturme 
für Eigengewicht und Wind. 

Die Abmessungen der Kegelschale mit dem Meridianschnitt Bild 14 sind nach einem 
vorhandenen Kühlturme festgelegt worden, der im vorliegenden Falle im Unterbau elastisch 
eingespannt sein würde. Die Schnittkräfte werden daher durch die Ergebnisse bei frei dreh- 
barer Lagerung und starrer Einspannung eingegrenzt. Die langwelligen Lösungsanteile 3, 4, 
die nach Zahlentafel 3 über die ganze Schale reichen, werden vernachlässigt. Die Störungen 


am unteren Rande lassen sich daher für 
und =27,8m 
unabhängig von den Bedingungen am oberen Rande berechnen. Die Integrationskonstanten 
C,—C, sind bei einem kräftefreien oberen Rande Null. 

1. Geometrische Grundlagen und Material. 
1„=126,68m; a=9%0', 0,00300. 


„618500, 1815 49878 - 


2. Die von der Belastung unabhängigen Zahlenwerte für n=0 bis n=3. 

a) Die kleinen Anteile 3, 4 der Lösungen für n=2 und n=3 erzeugen nach Bild 8 
periodische Abklingungsfunktionen. Die Abklingungszahlen für n=0, n=1 werden nach (35) 
und (37), die Abklingungszahlen für n=2, n=3 nach (16) streng berechnet und in Zahlen- 
tafel 2 mit den Ergebnissen des Näherungsverfahrens (31) verglichen. 


1/1944 
= ) 
E = 2'100'000 u=1/6. 
k = 15 = 7,5000 - 10 1393333. 
E 
1— 
1 
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Zahlentafel 2. Zahlentafel 3. 


Un bzw. r L, (m) 


58,827 | 58,864 | 1,5/0,5 ia 13,2 22,5 
59,151 | 58541 | 1,5/0,5 nn 13,4 23,2 
genau | 60,136 | 57,590 | 1,478 | 0,923 13,7 25,9 
genähert | 60,143 | 57,586 | 1,410 | 0.808 245,0 27,9 
genau | 61,842 | 56,089 | 2,974 | 2,505 13,9 31.0 
genähert| 61,855 | 56.080 | 2966 | 2,496 181,2 22,2 


b) Die Wellenlänge L, und das Verhältnis », der Schwingungsweiten werden mit (46) 


und (47) in Zahlentafel 3 angeschrieben. 
c) Die Faktoren f,, fr+, (r=1,3) werden nach Zahlentafel 1 berechnet. 


3. Schnittkräfte für Eigengewicht mit y =2,4 t/m?. 
a) Membranzustand für die äußeren Kräfte. 

_ 

y’ 


1 


N 


y tan 


Bild 15 zeigt u. a. die Schaulinien der Membrankräfte. 
Die Randwerte von Ringdehnung und Meridianverdrehung sind: 
& = + 2,667 8,010 10=®, 
b) Biegezustand- durch Randbelastung. 
Als Randbedingungen werden bei frei drehbarer, aber unverschieblicher Lagerung 


&p=— Eyyembran; My=0, bei starrer Einspannung — w — W'yembran VOr- 
geschrieben. Damit entstehen die Integrationskonstanten C, der Zahlentafel 4. 


Zahlentafel 4. 


Lagerung gelenkig (II) eingespannt (III) 
Rand- = — 2,667 - 10 
bedingungen | 9 w' = — 8,010-10=® 


Ö — 5,561 10-5 — 5,567 - 10-5 
— 9,086 - 10? + 3,793 


Zahlentafel 1 liefert außerdem die Abklingungsfunktionen der Komponenten des Spannungs- 
zustandes. 

ce) Überlagerung von Membran- und Biegezustand. 

Mit den Schnittkräften des Bildes 14 erfüllt die Schale die Auflagerbedingungen &, —0, 
M,=0 bei gelenkigem Anschluß und die Auflagerbedingungen &,=0, w’—=0 bei starrer 
Einspannung. 

4. Schnittkräfte unter dem Strömungswiderstand der Schale bei Wind mit dem Stau- 
druck pw» = 110 kg/m’. 

Die Druckverteilung p, über den Breitenschnitt ist nach Bild 16 in Anlehnung an die 
Göttinger Versuche über den Strömungswiderstand von zylindrischen Baukörpern gewählt 


worden: 
P: (9) = [— 0,50 + 0,50 cos # + 0,85 c0os29 +0,15 c0s3 9] : px. 


Die Druckverteilung p:, über den Meridian ist mit pz, (a) = pw c0S’a konstant. Äußere Kräfte 
sind demnach: 


(46) 


18 
T- 
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Bild 14. 


Bild 15. Schnittkräfte für Eigengewicht. 


2, 0,50 +0,50 c0s # +0,85 cos 29+0,15005839] -pucos’a, 


n=3 
P:= Z,cosnd; Zu =, Pw C08? a. 
n=0 
a) Membranzustand für die äußeren Lasten. 
Die inneren Kräfte sind mit den Bedingungen N,„= Nyyn=0 am oberen Rande y=|,: 


n=3 n=3 
Ny= Nynwsnd; Ny= Nyncosnd; Nyynsinn®. 
n—=0 n=1 
N —3sin’ a) — 3(n’— sin? + 


Bild 17 enthält die Schaulinien der summierten Membrankräfte.. Am Unterrande entsteht 
für die Harmonische » die Ringdehnung 


Ösina cosa sin?’ a+ un?) ulm’—sin’a) 7;+24 n? 


und mit v=v=v’=0 die Meridiandrehung 


‚Mech. 
| | 
| | 
N | | | 
- 80%’ | — Membranzustand 
= 
Nyn=—yZutana, 


— 
m 
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> 
97 
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— 

130 — 
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Bild 17. Schnittkräfte für Windbelastung. 


Bild 18 (rechts). Schnittkräfte für Windbelastung. Längen- 
maßstab in Grund- und Aufriß verschieden. 


Die Zahlenrechnung liefert für &y, b) Biegespannungen durch die Randkräfte. 
und mw, folgende Ergebnisse: Die Zahlenrechnungen werden getrennt für jede 
Harmonische n wie bei V 3b durchgeführt und dabei die 
Abklingungsfuuktionen ohne den kleinen Anteil 3,4 der 
Lösung gebildet. Um jedoch ihre Bedeutung für das Er- 
n=0| + 1,321-10=®| + 1,667- 10" gebnis zu untersuchen, werden die Integrationskonstanten 
» C,—C, und die Randwerte der Schnittkräfte und Ver- 
— 1,874-10°°| — 1,010-.10°%°  gchiebungen für n=2 und n=3 Zahlentafel 5 mit und 
2 z = ohne Berücksichtigung der kleinen Schwingung 3,4 an- 
EEE gegeben. Die Ergebnisse unterscheiden sich nach Zahlen- 
” - mr tafel 6 bis auf N, nur wenig und begründen damit die 
—3| — 1,891-10-#| — 2,327. 10-® 

Brauchbarkeit der Näherungslösung des Abschnittes IIID. 
Der Unterschied bei N, kennzeichnet die langwelligen Anteile der Abklingung als Längs- 
spannungszustände. Der Betrag dieser Randstörung ist jedoch so gering, daß er gegen die 
Membrankraft N, für die äußeren Lasten vollkommen zurücktritt. 


c) Überlagerung von Membran- und Biegezustand. 

Die Membran- und Biegeanteile der drei Schnittkräfte N,, Ny, Nyy werden zunächst 
für jedes » superponiert und damit die acht endgültigen Schnittkräfte jeder Belastungs- 
funktion n ermittelt. Die wirklich auftretenden Schnittkräfte werden nunmehr an jedem be- 
liebigen Schalenpunkte aus den einzelnen Gliedern der trigonometrischen Reihe gebildet. 
Die Bilder 17 und 18 enthalten die Schaulinien für einige Meridian- und Breitenschnitte. 
Die Bilder 19 und 20 zeigen den Einfluß der einzelnen Glieder der Fourierreihe auf M, 


und Myy- 


Un 


3 
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Darstellung von M, ın 
und Vertikalschnitten 


„7 


Bild 20. Drillungsmoment Myy = My [tem] infolge Windbelastung. Längenmaßstab in Grund- und Aufriß verschieden, 
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Zur praktischen Berechnung der Eigenwerte linearer 
Integralgleichungen. 


Herrn Friedrich Pfeiffer zum 60. Geburtstag am 2. Mai 1943 gewidmet. 
Von Friedrich Lösch in Rostock. 


Es wird ein Verfahren zur Berechnung der Eigenwerte linearer Integral- 
gleichungen dargestellt, das für den Fall des symmetrischen Kerns mit dem 
bekannten Grammelschen Verfahren übereinstimmt. Seine Brauchbarkeit 
auch bei unsymmetrischem Kern wird an dem Beispiel der Transversal- 
schwingungen eines einseitig eingespannten Stabs von veränderlichem 
Querschnitt gezeigt. 


Für die praktische Berechnung der Eigenwerte einer linearen Integralgleichung 


b 


mit symmetrischem Kern K(x.£) hat vor einiger Zeit Herr Grammel') ein Verfahren ent- 
wickelt, das sehr schnell zu besonders guten Näherungswerten führt. Im folgenden soll ge- 
zeigt werden, daß man zu den Grammelschen Gleichungen in einfacher Weise gelangt, 
indem man an die Hilbertsche Methode der unendlich vielen Variabeln?) anknüpft. Im 
Gegensatz zu dem von Herrn Grammel bei der Ableitung seines Verfahrens eingeschlagenen 
Weg läßt dieses Vorgehen nicht unmittelbar erkennen, daß die aus den Grammelschen 
Gleichungen entspringenden Näherungen durchweg obere Schranken der Eigenwerte darstellen, 
bietet aber andererseits den Vorteil, daß es ohne weiteres die Übertragbarkeit des Verfahrens 
auf Integralgleichungen mit unsymmetrischem Kern ergibt. Daß auch in diesem Falle das 
Verfahren gute Dienste leisten kann, wird an dem Beispiel der Transversalschwingungen eines 
am einen Ende eingespannten, am anderen Ende freien Stabes gezeigt. Man gelangt hier zu 
wesentlich besseren Resultaten, als sie z. B. das bekannte von Herrn Schwerin?) ausge- 
arbeitete Verfahren liefert. 


$ 1. Darstellung des Verfahrens. 


Wir legen den folgenden Betrachtungen die homogene Integralgleichung 
b 


a 
zugrunde, wobei wir über den Kern voraussetzen: 
(a) K (x, &) ist im Bereich a<x,2<b stetig; 
(b) K(x,&) ist im Bereich a<x,£<b gleichmäßig nach einem über 
orthogonalen und normierten Funktionensystem 


entwickelbar, es gilt also gleichmäßig 


b 
a 
Es liegt dann nahe, die Lösung der Integralgleichung durch einen Ansatz 


zu versuchen. Geht man damit in (1) ein, so wird man bei gliedweiser Ausführung der 
Integration — ob dieses Vorgehen gerechtfertigt ist, ist für das Weitere unwesentlich — zu- 
nächst auf die Gleichung 


1) R. Grammel: Ein neues Verfahren zur Lösung technischer Eigenwertprobleme. Ing.-Arch. Bd. 10 (1939), 
S. 35 bis 46. 

2) D. Hilbert: Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen. Leipzig und Berlin 1912, 
S. 174 bis 185. Vgl. auch G. Hamel: Integralgleichungen. Berlin 1937, S. 124 bis 130. 

3) E. Schwerin: Über Transversalschwingungen von Stäben veränderlichen Querscehnitts. Z. techn. Phys. Bd. 8 
(1927), S. 264 bis 271. 


3’ 


_ 
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= b 


geführt. Multipliziert man diese Gleichung nacheinander mit y, (x), y,(&), ..., yulz), »-. 
und integriert jeweils gliedweise über das Intervall a<x=<b, so ergibt sich das unendliche 


lineare Gleichungssystem 


wobei zur Abkürzung 

b b 

0 

a a 
gesetzt ist. Formal erscheint so die Frage nach der Existenz einer nicht identisch ver- 
schwindenden Lösung y(x) der Integralgleichung (1) zurückgeführt auf die Frage nach der 
Existenz einer nicht identisch verschwindenden Lösung (a,) des unendlichen Gleichungs- 
systems (5). 

Dieses Gleichungssystem, das der Hilbertschen Methode der unendlich vielen Variabeln 
zugrunde liegt‘), behandeln wir nun in der Weise, daß wir die endlichen Abschnitte des- 
selben betrachten und dann den Grenzübergang ausführen. Wir setzen demgemäß für 
n=1,2,3,.... das System der n homogenen linearen Gleichungen mit »n Unbekannten 


n 
an. Dieses endliche Gleichungssystem besitzt genau dann eine nichttriviale Lösung, wenn A 
Nullstelle der Determinante 


ist. Wir bezeichnen die n Nullstellen dieser Determinante mit 


wobei wir uns jede so oft aufgeführt denken, wie es ihrer Vielfachheit entspricht. 


Da bei unbegrenzt wachsendem n das Gleichungssystem (7) in das unendliche System (5) 
übergeht, so lassen unsere Überlegungen vermuten, daß dabei unter geeigneten Voraussetzungen 
die Werte (9) in die Eigenwerte der Integralgleichung (1), d. h. die Nullstellen der Fred- 
holmschen Determinante von (1) übergehen. Dies trifft in der Tat unter den zu Anfang 
gemachten Voraussetzungen über den Kern K’(x,£) zu. Es gilt nämlich der folgende 


1. Satz: Unter den Voraussetzungen (a) und (b) über den Kern K (x, &) gilt für die 
Determinante D, (4) gleichmäßig in jedem beschränkten Bereich der komplexen i-Ebene 


wo D (4) die Fredholmsche Determinante der Integralgleichung (1) bedeutet. Damit 
konvergieren für n>x die Nullstellen von D„(4) nach Lage und Ordnung gegen die 
Nullstellen von D (4), d. h. gegen die Eigenwerte der Integralgleichung (1). 


Beweis: Der Beweis dieses Satzes ergibt sich sehr einfach im Anschluß an bekannte 
Untersuchungen von E. Goursat und H. Lebesgue°) über die Approximation von Kernen 
durch ausgeartete Kerne. Wir entnehmen der Entwicklung (3) des Kerns K (x, &£) die Folge 


ausgearteter Kerne 
n 


4, Vgl. e.2) 
. 5) E. Goursat: Sur un cas @lementaire de l’equation de Fredholm. Bull. Soc. math. France Bd. 35 (1907), 
S. 163 bis 173. H. Lebesgue: Sur la möthode de M. Goursat pour la resolution de l’&quation de Fredholm. Bull. 
Soe. math. France Bd. 36 (1908), S.3 bis 19. Vgl. auch R. Courant-D. Hilbert: Methoden der mathematischen Physik 1. 


Berlin 1924, S. 124 bis 138. 
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die nach Voraussetzung (b) in a<x, £<b gleichmäßig gegen K (x, &) konvergiert. Setzt man 
die Lösung der zum n. ausgearteten Kern gehörigen Gleichung 
b 


a 


in der Form 


n 


an, so erhält man für die Konstanten 5{” in bekannter Weise das lineare Gleichungssystem 


n 


v=1 


wobei zur Abkürzung 


b b 


gesetzt ist. Die eingangs genannten Untersuchungen von E. Goursat und H. Lebesgue 
besagen dann, daß für n>x die Determinante des Gleichungssystems (14) gegen die Fred- 
holmsche Determinante D (4) der Integralgleichung (1) konvergiert, und zwar gleichmäßig in 
jedem beschränkten Bereich der komplexen 4- Ebene. 

Weiter stimmt nun aber das Gleichungssystem (14) mit dem System (7) überein, denn 
den Darstellungen (3) und (11) der Kerne K(x, &) bzw. K‘'” (x, &) entnimmt man, indem man 
nur die Orthogonalität der Funktionen y, berücksichtigt, daß 


K)=K,, für 


gilt. Damit ist das Bestehen der Limesrelation (10) schon bewiesen. 

Auf Grund dieser Beziehung (10) vollzieht sich nun der Nachweis, daß die Nullstellen Pr 
von D„(4) für n>x nach Lage und Ordnung gegen die Nullstellen von D (A) konvergieren, 
durch übliche Schlüsse: Man hat nur zu beachten, daß nach dem Satz von Rouche die 


beiden Funktionen 


D(ü) und Di) 


im Innern eines beliebigen Kreises der komplexen i-Ebene, auf dessen Peripherie D (4) =F0 
ist, gleich viele Nullstellen aufweisen, sobald längs der Peripherie 


ID, 
was nach (10) für genügend große n gewiß der Fall ist. 


Der im vorstehenden angewandte Grenzprozeß liefert, wie weiter bemerkt sei, nicht 
nur die Eigenwerte, sondern auch die Eigenfunktionen der Integralgleichung (1). Der Kürze 
halber soll dies hier nur unter der Annahme, daß die Eigenwerte der Integralgleichung sänıt- 
lich einfach sind, auseinandergesetzt werden. Im allgemeinen Falle hat man ganz analog 


zu schließen. 
Wir greifen eine (auf Grund unserer Annahme einfache) Nullstelle A, der Fredholm- 


schen Determinante D (4) heraus. Sodann bezeichne 

eine gegen A, konvergierende Folge von Nullstellen der Determinanten 
D, (A), D, (A), D, (A), 


Für genügend großes n sind auch diese Nullstellen einfach. Es gehört dann zu A jeweils 
eine Lösung der Gl. (7): 


( ) n) 
ah), ... a‘) 


die überdies eindeutig bestimmt ist, wenn wir sie gemäß der Bedingung 


n 


fech. 
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normiert denken. Mit dieser Lösung bilden wir — entsprechend dem Ansatz (4) — die über 
a<x<b normierte Funktion 


n 


Und nun gilt der 


2, Satz: Der Kern K (x,£) erfülle die Voraussetzungen (a) und (b) und besitze über- 
dies nur einfache Eigenwerte. Ist dann 4, irgendein Eigenwert, so yilt für die durch (16) 
erklärte Funktion y” (x) gleichmäßig im Intervall azx=b 


im (@)=y(), 
n>x» 


wo y(x) die zum Eigenwert }, gehörige normierte Eigenfunktion des Kerns K (x, £) bedeutet. 


Beweis: Wie aus dem Beweis des 1. Satzes bekannt ist, stimmt das Gleichungs- 
system (7) mit dem System (14) überein. Der Vergleich von (16) mit dem Ansatz (13) zeigt 
daher, daß die durch (16) gegebene Funktion y'” (x) die zum Eigenwert A" gehörige nor- 
mierte Eigenfunktion des ausgearteten Kerns K"' (x, £) darstellt. Auf Grund der genannten Unter- 
suchungen von E. Goursat und H. Lebesgue folgt dann die Behauptung unmittelbar aus der 
Tatsache, daß die ausgearteten Kerne den Kern K (x,£) gleichmäßig approximieren. 

Das aus den beiden vorstehenden Sätzen entspringende Verfahren zur Berechnung der 
Eigenwerte und Eigenfunktionen deckt sich für symmetrische Kerne im wesentlichen mit dem 
von Herrn Grammel angegebenen Verfahren‘). In der Tat stimmen die Grammelschen 
Gleichungen, falls nur wie hier die Orthogonalität und Normiertheit der Koordinatenfunktionen 
y, angenommen wird, genau mit den obigen Gl. (7) überein. Herr Grammel gewinnt die 
Gleichungen unter etwas anderen Voraussetzungen aus einem Variationsprinzip und hat damit 
gegenüber der vorliegenden Behandlung den großen Vorteil, zu wissen, daß die aus seinen 
Gleichungen fließenden Näherungswerte durchweg obere Schranken der Eigenwerte darstellen. 
Es ist jedoch u. U. zweckmäßig, das Verfahren auch für unsymmetrische Kerne anzuwenden, 
und hier versagt der Weg über das Variationsprinzip. In der Tat zeigt die in $ 2 gegebene 
Anwendung, daß im Falle eines unsymmetrischen Kerns die mit Hilfe der Gl. (7) berechneten 
Näherungen der Eigenwerte nicht notwendig obere Schranken derselben darstellen. Bei sach- 
gemäßer Anwendung des Verfahrens sind aber auch in diesem Falle die sich ergebenden 
Näherungswerte sehr gut, was ja für symmetrische Kerne durchaus bekannt ist’). 


$ 2. Anwendung des Verfahrens. 
Transversalschwingungen von Stäben veränderlichen Querschnitts. 


Wir gehen nun noch etwas näher auf die praktische Anwendung des aus den Betrach- 
tungen von $ 1 entspringenden Verfahrens zur Berechnung der Eigenwerte einer linearen 
Integralgleichung (1) ein. Hier kommt es darauf an, daß schon die ersten Schritte gute 
Näherungen liefern. Unter diesem Gesichtspunkt sind die Funktionen y, zu wählen. Nehmen 
wir nun an, es seien 


die Eigenwerte der zu behandelnden Integralgleichung und 
stelle das zugehörige normierte, orthogonale System von Eigenfunktionen dar, so gilt durchweg 
\ Kia, =, y,(@), 
a 


und damit reduziert sich das Gleichungssystem (5) auf 


Durch Nullsetzen der zugehörigen Determinanten 
6) Vgl. 1. c.t). 


?) Vgl.l.e.t) sowie L. Collatz: Genäherte Berechnung von Eigenwerten. Z. angew. Math. Mech. Bd. 19 (1939), 
S. 224 bis 249 und S. 297 bis 318. 
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erhält man daher nacheinander exakt die Werte 


Diese Überlegung läßt annehmen, daß man für die Eigenwerte gute Näherungen er- 
halten wird, wenn man für die ersten Funktionen y, einigermaßen gut geschätzte Näherungen 
der entsprechenden Eigenfunktionen annimmt. Für symmetrische Kerne kommt man so zu 
dem beim Grammelschen Verfahren üblichen Vorgehen, das in der Tat sehr gute Werte 
liefert. Im folgenden soll an einem Beispiel gezeigt werden, daß man auch bei unsymmetri- 
schem Kern günstige Resultate erwarten kann. ö 


Als Beispiel wählen wir das Problem der Transversalschwingungen von Stäben ver- 
änderlichen Querschnitts. Dieses Problem ist unter Heranziehung der Theorie der Integral- 
gleichungen bereits von Herrn Schwerin*) ausführlich behandelt worden. Seine Methode be- 
steht darin, den Kern der zuständigen Integralgleichung durch die ersten Glieder seiner 
Fourier-Entwicklung zu approximieren und dann die Theorie der Integralgleichungen mit 
ausgeartetem Kern anzuwenden. Wir werden im folgenden das Problem nach der oben ent- 
wickelten Methode behandeln und uns dabei genau derselben Annahmen wie Herr Schwerin 
bedienen. 

Es bezeichne für einen am einen Ende («=0) eingespannten, am andern Ende («=1) 
freien Stab: 

F(x) die Querschnittsfläche in x, 
F, die Querschnittsfläche in 2=0, 
J(x) das Trägheitsmoment in «, 
J, das Trägheitsmoment in x=0, 
E den (konstanten) Elastizitätsmodul, 
o die (konstante) Dichte des Stabmaterials. 


Sowohl für das reziproke Trägheitsmoment als auch für den Stabquerschnitt werde an- 
genommen, daß sich ihr Verlauf hinreichend genau durch eine Ausgleichsparabel 
1 
J (x) 
darstellen lasse. Dann besteht für die Eigenfrequenzen p und die zugehörigen Schwingungs- 
formen y(x) des Stabs die Integralgleichung 
1 


0 
mit 
x|5+ 6 + 12 für 
G(&,d)= 2 3 4 3 5 (18a) 
und 
8b). 


Wie man weiß, ist der Kern der Integralgleichung (17) positiv definit; die Eigenwerte 
denken wir nach wachsender Größe numeriert: 


Um nach der in $ 1 entwickelten Methode*®) Näherungen für die niedersten Eigenwerte zu 
bekommen, wählen wir für y, (x) ein normiertes Polynom möglichst niedriger Ordnung, das 
die für die Schwingungsformen (x) geltenden Randbedingungen 

erfüllt. Dies leistet das Polynom 4. Grades 

it x* 
mi 
a,=3,946761, 23631174, y,= 0,6577. 


s) Vgl. 1. e.3). 
®) Da der Kern symmetrisierbar ist, so läßt sich direkt das Grammelsche Verfahren anwenden. Es empfiehlt 
sich jedoch, an die unsymmetrische Form der Integralgleichung anzuknüpfen, damit sich die Elemente Ku» der Deter- 


minante D„(4) bequem in den Parametern a;, «aa, ausdrücken. 
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Für y,(x) wählen wir sodann das normierte Polynom 5. Grades, das die Randbedingungen (19) 
erfüllt und zu y, (x) orthogonal ist: 


a,—=32,162440, P,=— 81,294019, y,=65,2127%9, d,= — 17,%5718. 


Entsprechend wählen wir schließlich für y, (x) das normierte Polynom 6. Grades, das die 
Randbedingungen (19) erfüllt und zu y, (x) und y,(x) orthogonal ist: 


mit 


mit 
a, = 117,604528, — 568,165629, y,=%9,328359, — 686,886776, = 179,985998. 
Brechen wir an dieser Stelle ab, so sind wir in der Lage, die Determinante 
D, (4) 

aufzustellen, womit zugleich D,(4) und D,(A) bekannt sind. Berechnen wir ihre Elemente 
Ku, nach den Formeln (6), so erhalten wir dieselben in der Form 

(20), 


13) »(19) (33) 
+ 


wobei die Größen Ko) die in der nachfolgenden Zahlentafel angegebenen Werte besitzen. 


Zahlentafel zur Berechnung der Größen Ku». 


„| 102 | 10: KU) | 10: | 10: | 10° | 10° 10: 10: 102 
| || | | | | | | 


80874 1,5590 0,4943 | 6,5051 5,4305 | 1,2848 | 1,0904 0,4154 | 0,3575 
0,1115 | — 0,2191 | — 0,1186 — 0,1218 0,1205 | — 0,1898 | — 0,1668 | — 0,1030 | — 0,0907 


0,0014 0,0076 0,0180 | 0,0042 0,0062 0,0099 0,0108 0.0170 0,0160 


0,1115 0.2191 0,1186 | — 1,3143 | — 1,8034 | — 0,3665 0,4249 | — 0,1470 | — 0,157 
02061 | 0,0865 0.0474 0,168 0,1065 | 0,0868 0,0853 0,0482 | 0,0473 
0,0158 | — 0,0205 | — 0,0150 — 0,0145 — 0,0124 — 0,0154 | — 0,0130 — 0,0131 | — 0,0119 
0,0014 0,0076 | 0,0180 0.2590 0,5928 | 0,0920 0,1601 0,0466 0,0660 
0,0158 | — 0,0205 | — 0,0150 | — 0,0509 | — 0,0576 | — 0,0330 | — 0,0414 | — 0,0192 | — 0,0259 


0,0270 0,0127 | 0,0095 | 0,0167 | 0,0141 | 0,0133 | 0,0095 | 0,0087 0,0075 


Um einen Anhaltspunkt für die Güte der damit erreichten Näherung zu gewinnen, be- 
handeln wir zunächst den Stab von konstantem Querschnitt, für den die genauen Werte 


i,=12,36%4, i,=485,66, A, 3806 


(jeweils auf die angeschriebenen Stellen abgerundet) bekannt sind. In diesem Falle wird 
einfach: 


2 
2 
2 


Ku 
und man erhält nacheinander 


aus D,(A)=0 die Näherung: 4,» 12,3648; 
aus D,(A)=0 die Näherungen: A, 12,3624; A, 488,70; 
aus D,(4)=0 die Näherungen: A, » 12,3624; 485,50; A, 3840. 


Demgegenüber erhält man bei demselben Arbeitsaufwand nach dem Schwerinschen Ver- 
fahren, wenn die dort zur Verwendung kommenden, den D, (4), D,(4), D, (A) entsprechenden 


Determinanten mit 
4, (A), A,(4), A, (4) 


bezeichnet werden, nacheinander '°) 


aus 1,(4)=0 die Näherung: 4, 
aus 4,(4)=0 die Näherungen: 4, = 12,71; 4,» 691,0; 
aus 4,(4) = 0 die Näherungen: 4, 12,355; A, 520,33; A, 6917. 

10) Bei der Ermittlung dieser Näherungen sind bereits einige Versehen berücksichtigt, die sich in dieSchwe- 
rinsche Arbeit eingeschlichen haben. Wie ich einer freundlichen Mitteilung von Herrn Pfeiffer entnehme, muß es dort 
in der Tabelle zur Berechnung der Größen Ayın neißen: 

Aoas = + 6,971 statt + 6,911; G33 = + 0,201 statt + 2,201; La3 = + 0,399 statt + 0,296. 
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Man erkennt, daß selbst die aus A, (A)=0 errechneten Näherungen für A, und 4, noch erheblich 
schlechter sind, als die oben für A, und 4, aus D,(A)=0 bzw. D,(4)=0 erhaltenen Werte. 


Als zweites Beispiel werde der durch 


> =1 + 0,79646 x + 0,99452 0,57197 x — 0,02357 x? 


charakterisierte Stab behandelt. Für ihn ist 
a,=0,79646, a,=0,%94532, a,=—0571%7, a, =— 0,02357; 
mit diesen Werten sind die Größen K., nach (20) an Hand der Zahlentafel zu berechnen. 


Man erhält dann nacheinander 


aus D,(A)=0 die Näherung: 


19,53; 


aus D,(A)=0 die Näherungen: A, = 19,6181; A, » 510,75; 
aus D,(A)=0 die Näherungen: 4, » 19,6184; A, » 511,80; A, » 3779. 


Das Schwerinsche Verfahren liefert dagegen mit demselben Arbeitsaufwand 


aus 4, (A)=0 die Näherung: 


A, = 11,64; 
aus A,(4)=0 die Näherungen: A, » 19,97; 


= 69,2; 


aus A,(4)=0 die Näherungen: 4, » 19,629; A, » 540,70; A, = 6649. 
Das Bild ist also genau dasselbe wie im Falle des Stabes mit konstantem Querschnitt. 520 


Eingegangen am 18. 10. 1943. 
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Zum Determinantenverfahren von Chiöd. 
Das Verfahren von Chiö zur zahlenmäßigen Be- 
rechnung von Determinanten, empfohlen in E. T. 
Whittaker und G. Robinson, The calculus of ob- 
servations, 1. Aufl. London 1924, 2. Aufl. London 
1926, Kap. V, besteht darin, die gegebene Deter- 
minante zu reduzieren auf das mit 
Vorzeichen (—1)P+@ versehene Produkt aus 
einem beliebig (wenn vorhanden, gleich 1) ge- 
wählten „Leitelement“ «»g mit einer Determinante 
nächstniedriger Ordnung, deren Elemente folgen- 
dermaßen zustandekommen: In der ursprünglichen 
Determinante streicht man die pte Zeile und die 
4te Spalte des Leitelements  apg und fällt von 
jedem übrigbleibenden Element a;, Lote auf die 
Streichungslinien. Dann nimmt man das Produkt 
der beiden getroffenen „Fußpunktelemente“ apk 
und a@ig, teilt es durch das Leitelement a», und 
zieht den so gefundenen Ausdruck a. 


pq 
dem alten Element @;, ab. Das neue Element heißt 
also 


(N). 
Apaq 
Diese Reduktionsregel für Determinanten ist 
weiter nichts als eine auf Determinanten zu- 
geschnittene Fassung des Gaußschen Eliminations- 
verfahrens für lineare Gleichungssysteme 


n 

Denn nach Gauß eliminiert man bekanntlich eine 
Unbekannte, z. B. die qte Unbekannte x. dadurch, 
daß man eine Gleichung, etwa die pte Gleichung 

DL Apk (3), 

k=1 


als „Eliminationsgleichung“ der Reihe nach mit 
geeigneten Faktoren so multipliziert, daß bei 
Addition zu jeder der übrigen Gleichungen zy her- 
ausfällt. In (2) hat zz die Vorzahl aig, in (3) die 


ai 
Vorzahl a»,. Also muß (3) mit _—_ multipliziert 

werden. Dabei entsteht das neue, die Unbekannte 
2q nicht mehr enthaltende Gleichungssystem 


n 
a; 


dessen Koeffizienten offenkundig der Bildungs- 
vorschrift (1) genügen. 

In der Determinante von (2) werden bei der 
Gaußschen Elimination lediglich Zeilen kombi- 
niert und dadurch in der qten Spalte Nullen her- 
beigeführt bis auf apg. Die Determinante von (2) 
wird also gleich (—1)P+9 a); mal der Deter- 
minante von (4), was mit der Regel von Chiö 
gleichbedeutend ist. 


Darmstadt. A. Walther. 499 


Über die Dämpfung von Schwingungen 
bei zeitlich veränderlichen Kräften. Ein 
System gekoppelter Freiheitsgrade %,,..., Ym möge 
einem System von m linearen Differentialgleichun- 
gen mit der Zeit £ als unabhängiger Variabler 
genügen. Sind die Koeffizienten zeitlich konstant, 
so entscheiden allein die Realteile der Wurzeln der 
in bekannter Weise gebildeten charakteristischen 
Gleichung des Systems über die Dämpfung der 
Schwingungen, die die 4. ausführen. Welche Ab- 
änderung der Regel tritt ein, wenn die Koeffizienten 
des Systems mit 2 langsam veränderlich sind ? 

Zur Beantwortung dieser Frage bringe man das 
System in die Gestalt 


n 


— 


mit nZ> m, wozu i. a. nötig sein wird, einige Ab- 


leitungen der Freiheitsgrade %,....,/„ nach der 
Zeit als neue Variable Ya+ı,-..,4n einzuführen. 
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Das System (1) läßt sich unter Verwendung der 
Matrizensymbole 


in der Gestalt 


schreiben. Führt man an Stelle von y ein System 3 
neuer unabhängiger Funktionen 2,,...,2„ mittels 
der Gleichungen 


n 


bzw. wieder in Matrizenschreibweise 


- 
06), 


ein, so gilt 
mit 


wobei 
d, 
uv = dt uv 
bedeutet. 

Nach Erwähnung dieses bekannten Sachverhaltes 
wähle man für R insbesondere eine solche Matrix, 
daß 


In 


eine Diagonalmatrix ist. Die /, sind die Wurzeln 
der charakteristischen Gleichungen der Matrix N, 
wir setzen sie hier als von einander verschieden 
voraus, damit (7) überhaupt möglich ist. Die 
Matrix AR kann hierbei noch mit einer beliebigen 
Diagonalmatrix ® links multipliziert, d. h. durch 
DR ersetzt werden. Dadurch kann erreicht 
werden, daß eine Zeile von R!, etwa die erste, 
zeitlich konstante Koeffizienten besitz. Wegen 

ist nun die Dämpfung der Bewegung von %, (£) die- 
selbe wie der Bewegung der z,(£); es genügt also, 
letztere zu untersuchen. 

Die z, bestimmen sich aus dem Differential- 
gleichungssystem (5), dessen Matrix nach (6) und (7) 


d= +6, €=RR... 
An 

ist, wobei noch anzumerken ist, daß die Koeffizi- 

enten der Matrix € unter anderem von den 


Fans homogen und linear abhängen; sie 


verschwinden also, wenn letztere verschwinden, 
und sind klein, wenn letztere klein sind. 

Sind diese Größen sehr klein, so erhält man bei 
ihrer Vernachlässigung CE=0 und die Lösungen 


dur 


(9) 


2,e) 


von (5), dabei bedeuten c, Integrationskonstanten. 
Bei langsam veränderlichen Koeffizienten au» ist 
(9) die erste Näherung an die Lösung. 

Um die zweite Näherung zu erhalten, trans- 
formiere man nochmals auf Diagonalform: 


(10) 


und führe neue Veränderliche 2,,...,2„ ein, welche, 
in einer einspaltigen Matrix 3 angeordnet, den 
Gleichungen a 


und 
(12) 
genügen, wobei 


+T, (13) 


ist. Die Elemente der Matrix € bestehen nun, wie 
aus dem folgenden hervorgeht, aus linearen Kom- 
binationen der zweiten Ableitungen ä.» und der 
Produkte üauräaww und sind folglich von höherer 
Ordnung klein als die Koeffizienten von €. Die 
durch Vernachlässigung von € entstehende Lösung 


e! 


(14) 


von (5) ist demnach besser als (9); es ist die ge- 
suchte zweite Näherung. 

Zur Auswertung von (14) braucht man die Werte 
der A, sowie die Matrix ©. Letztere wird sich 
von der Einheitsmatrix E nicht wesentlich unter- 
scheiden, wenn die Elemente von E hinreichend 
klein sind, wie angenommen werden sollte. Man 
setze also 


S=-E+T, 
dann ist nach (8) und (10) 
L+E+TL— 
also gilt für die Elemente Zu» von T 
t,=0, tur(l,— für 


und 


Die Formeln (11) und (14) ergeben dann 


e) Au + dt 


Au — 


(15). 


Im Sinne unseres bisherigen Verfahrens, 
Äur—=üuräww®—=( zu setzen, haben die Koeffizi- 
enten dieses Systems als Konstanten zu gelten. 
Folglich sind die Funktionen 2,(£) in gleicher 
Weise gedämpft wie die Funktionen (14). 

Das hiermit gewonnene Ergebnis kann folgender- 
maßen zusammengefaßt werden: 


Regel: Es sei das System (1) linearer ho- 
mogener Differentialgleichungen ersten 
Grades vorgelegt, dessen Koeffizienten 
von £ schwach abhängen. Dann ist die 
i-te (z.B. 1-te) gesuchte Funktion %; (z. B. %,) 
in erster bzw. zweiter Näherung eine 
lineare Kombination der Funktionen 


dt 
bzw. ; 


dabei bedeuten A, die Wurzeln der cha- 
rakteristischen Gleichung der Matrix 
A=(au») und c,, die Diagonalenelemente 
der Matrix 
d 

wobei R der Gl. (7) genügt und dabei so 
normiert ist, daß die Elemente der i-ten 


:he, 
den 


11) 


13) 
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(z. B. I-ten) Zeile von R”! von £ unab- 
hängig sind. 


Beispiel: Eine Differentialgleichung 
zweiter Ordnung läßt sich als System so schreiben: 
Ya, 
»=—-by— ap. 


Es ist nun 


22 


Die Bewegung von %, verläuft also gedämpft, wenn 
die Realteile von 


a 
negativ sind; die Beträge dieser Realteile sind die 


logarithmischen Dekremente der Bewegung. 
Handelt es sich um echte Schwingungen: 


und o und w reell, so sind diese stabil, sobald 


92 


ist, anderenfalls labil. 


Göttingen. Martin Eichler. 510 


Uber einen mit Nilfe seines Schattens 
geienkten Integraphen'). Bild 1 zeigt zwei 
windschiefe Gerade Z, G in Auf- und Grundriß. 
Z ist lotrecht, G unter 45° geneigt. bd ist das 
Gemeinlot der beiden und habe die Länge /; c sei 
der Punkt auf G, dessen Aufriß c” auf Z” liegt. 
Dann ist 

wo a der Winkel des Gemeinlotes mit der Rich- 
tung des Aufrisses ist. Wir werden meist !=1 
nehmen. 


14 


Bild 1 und 2. 


1) Vortrag auf der Kriegstagungder Deutschen Mathema- 
tiker-Vereinigung 1943 in Würzburg. 


Nun bringen wir G und Z als gespannte Fäden 
an einem zweirädrigen Karren so an, daß die 
Achse seiner Räder in der lotrechten Ebene durch 
bd liegt. Die Ebene, auf der er fährt, sei die 
Grundrißebene //;; ferner schreibe er mit einem 
federnd auf //, drückenden Stift im Schnittpunkt 
von Z mit I/II. Seine Deichsel liege in der durch 
Z gehenden zu G parallelen Ebene und habe an 
ihrer Spitze einen auf //, gleitenden Fuß. Er- 
zeugen wir dann noch den Aufriß des Geraden- 
paares als Schatten in einer Parallelbeleuchtung, 
so haben wir einen Grundintegraphen, d. h. ein 
Gerät zum Zeichnen von Integralen (f(z)de. 
Um nämlich ein solches Integral zu zeichnen, 
bringen wir die Kurve z=f(z) so in die Aufriß- 
ebene, daß d’”’ und b” in der etwa nach rechts ge- 
ichteten x-Achse liegen. In der Grundrißebene 
verwenden wir dasselbe z, während die %-Achse 
normal zur Aufrißebene nach hinten gerichtet sei. 
d’ habe in I//, die Koordinaten =, y. Wir ziehen 
den Karren an seiner Deichsel so dahin, daß der 
Schnitt der Schatten von G und Z auf der Kurve 
(x) wandert. Dann ist beständig y’=1f 
also 


y=lii(e)de. 


Die einzige Schwierigkeit dabei ist, daß es an- 
scheinend noch keine hierfür brauchbare Parallel- 
beleuchtung gibt. Die Sonne liefert sie jedenfalls 
nicht ohne weiteres. Denn sie hat einen schein- 
baren Durchmesser von %°, d. h. die Sonnen- 
strahlen weichen von einer mittleren Richtung 
nach allen Seiten um Winkel bis zu %° ab. 

Wir können Z” und G” aber auch ohne Parallel- 
beleuchtung auf folgende Weise als Schatten be- 
kommen. Wir bringen auf dem Karren zwei elek- 
trische Glühbirnen, in denen je ein gerader Faden 
leuchtet. so an. daß die leuchtenden Fäden in den 
Verlängerungen der Fäden Z und G liegen. Diese 
T.ampen blenden wir so ab, daß sie die Aufriß- 
ebene nicht unmittelbar, sondern durch einen zur 
Aufrißebene parallelen ebenen Spiegel beleuchten. 
Die Spiegelbilder G* und Z* werfen dann von G 
und Z auf //,» xerade die Schatten G” und Z”. 
Zwar ist G” von Z*, Z” von G* beleuchtet, sie 
sind aber trotzdem sehr deutlich, vorausgesetzt. 
daß sonst kein Licht auf 77, fällt, weil das Auge 
einen Sprung in der Helligkeit von 50% auf 
100% sehr gut wahrnimmt. Natürlich wirft Z* 
auch einen Schatten von G, G* einen Schatten 
von Z: die sind aber so blaß, daß sie nicht stören. 


Bringt man das Geradenpaar Z, G so auf dem 
Karren an. daß die Achse seiner Räder mit db 


den Winkel +5 einschließt, dann wird 


Ist und /=1. 


wo k=tgß gesetzt ist. 


dann wird 


1 
Der Integraph zeichnet dann also das Kehrwert- 
integral 
Statt der Geraden Z und G kann man auch G 
und die aus G@ durch Drehung um Z um 180° ent- 


stehende Gerade H verwenden. Das ist bei kleinem 
! zweckmäßig, weil es die Genauigkeit vergrößert. 
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Bei größerem ! bekommt der Karren dadurch aber 
eine unhandliche Breite. Es ist unwesentlich, daß 
G (und damit H) unter 45° geneigt ist. Eine 
Änderung der Neigung wirkt wie eine Änderung 
von 1. 


Man kann diesen Integraphen auf zwei Arten 
leicht zu einem Momenten- und Stieltjesintegraphen 
verallgemeinern. Verfügt man nämlich über eine 
hinreichend saubere Parallelbeleuchtung, dann 
kann man G durch eine in einer lotrechten Ebene 
& liegende Kurve K ersetzen (Bild 2). « und z 
seien Koordinaten in ge, wobei 2 wie oben die 
Höhe über einem Punkt d von Z bedeutet. Das 
Lot von d auf e treffe e im Punkt b. u sei der 
Abstand des Punktes (w,z) von der Parallelen zu 
Z durch 5, und zwar habe ein in der zu //» nor- 
malen Ebene durch Z liegender Punkt c posi- 
tives v4, wenn a positiv ist. Der Einfachheit halber 
sei Z=1 und die Radachse zu e normal. Wir 
ziehen den Karren so an der Deichsel dahin, daß 
der Schnittpunkt der Schatten von Z und K auf 
der Kurve z=f(x) wandert. Hat K ine die 
Gleichung g (z2), dann wird 


d. h. 


(2). 


Unter diese Form fallen die Momentintegrale 
Schreibt man (2) in der Form 
y=\g(z)de, wo ist, 


so hat man ein Stieltjesintegral. Diese Integra- 
tıon kann im Gegensatz zu anderen Stieltjesinte- 
graphen ein Arbeiter ausführen ?). 

Bei Verzicht auf eine Parallelbeleuchtung hat 
man folgende Möglichkeit. Man verwendet als 
schattenwerfende Kurven wieder Z und G, als 
Lichtquellen wie oben Z* und G*, stellt aber mit 
Hilfe eines Funktionsgetriebes einen funktionellen 
/usammenhang zwischen den obengenannten 
Winkeln a und her. Er sei in der Form a + ß 
— h(a) geschrieben. Dann wird 


=tg 
wo tgh (arctg 2) = 9 (2) gesetzt ist. 


Innsbruck. L. Vietoris. 5231 


2) Vgl. etwa Mever zur Capellen, Mathematische 
Instrumente. Leipzig 1941, S. 188. 
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Mitteilungen aus dem Institut für 
Aerodynamik der Eidgen. Techn. 
Hochschule Zürich. Nr. 8. Herausgegeben 
von Prof. Dr. J. Ackeret. 72 S. Zürich 1943, 
Verlag A. G. Gebr. Leemann & Co. 


Das Heft beinhaltet vier Aufsätze: (1) Abdruck 
einer schon früher veröffentlichten Beschreibung der 
Windkanäle des Züricher Institutes (J. Ackeret); 
(2) Beschreibung eines Gerätes zur Messung tur- 
bulenter Schwankungen in Windkanälen mit Hilfe 
der Hitzdrahtmethode, das gegenüber dem be- 
kannten Gerät von Dryden vereinfacht ist (G. 
Dätwyler); (3) Einzelheiten der Schlierenoptik 
des Züricher Überschallkanales (P, de Haller). 
(4) An einem Tragflügel ergab bei der. Widerstands- 
messung ein Vergleich der Impulsmethode mit der 
Wägung recht gute Übereinstimmung, wenigstens 
für kleine Auftriebswerte. Bemerkenswert ist die 
erstmals mitgeteilte Bestimmung des laminar- 
turbulenten Umschlagspunktes mit Hilfe eines. an 
ps Staurohr angeschlossenen Hörrohres, wie es 

rzte verwenden; das Einsetzen der Turbulenz wird 
durch charakteristische Geräusche (Knattern) an- 
gezeigt (W. Pfenniger). 


Göttingen. H. Schuh. 514 


Dr.-Ing. Dr. techn. E. RAUSCH, a. pl. Prof. a. d. 
Techn. Hochschule Berlin, Maschinenfun- 
damente und andere dynamische 
Bauaufgaben,3. Teil: Rahmenfunda- 
mente bei hoher Maschinendrehzahl 
(insbesondere Dampfturbinenfunda- 
mente), Aufstellung von Maschinen 
in Gebäuden, Werkzeugmaschinen, 
Schäden an Maschinenfundamenten 
undErschütterungsschädenanderer 
Art, andere dynamische Aufgaben 
im Bauwesen nebst Ergänzung zum 1. und 
2. Teil des Buches. 365 S. m. 377 Bildern. Berlin 
1942, VDI-Verlag. Preis brosch. 35 M. 

Der Verfasser, der in dem 1. Teil des Buches 
ragen Gesichtspunkte für die bauliche Durch- 
bildung von Maschinenfundamenten und im 2. Teil 


Ausführungsbeispiele für Block- und Kastenfun- 
damente bei hin- und hergehenden oder um- 
laufenden Maschinen angegeben hat, behandelt 
in diesem 3. Teil die Rahmenfundamente für 
Maschinen mit hoher Drehzahl und beschließt die 
Arbeit mit wertvollen Hinweisen über die Auf- 
stellung von Maschinen, über Bauschäden an Ma- 
schinenfundamenten, Erschütterungsschäden und 
andere Bauaufgaben, deren Lösung zu dynami- 
schen Untersuchungen führt. Das Buch ist für 
den Bauingenieur geschrieben, der Belehrung über 
die Anwendung der einfachen dynamischen Be- 
rechnungen sucht, mit denen die elasto-kineti- 
schen Eigenschaften der Maschinenfundamente be- 
urteilt und als Grundlage für die Bemessung und 
bauliche Durchbildung der Baukörper verwendet 
werden. Er schöpft hierbei aus der eigenen, 
reichen Erfahrung und berichtet mit ausführlichen 
Zahlenrechnungen und Zeichnungen vor allem 
über die allgemeine’ Anordnung, die statische Be- 
rechnung, die Bemessung der Rundstahlbewehrung 
und über die Bauausführung von Dampfturbinen- 
fundamenten. Die eigentlichen dynamischen Ge- 
sichtspunkte treten dabei zurück, obwohl man in 
diesem Zusammenhang gern die Bestätigung der 
Rechnung durch Einmessen der Resonanzlagen 
beim Hochfahren der Maschinen durch Amplitu- 
denkurven erhalten hätte. 

Hieran schließen sich Angaben über die Auf- 
stellung von Maschinen in Gebäuden und über 
Bauschäden an Fundamenten, die durch dyna- 
mische Wirkungen in Verbindung mit mangel- 
hafter Bauausführung hervorgerufen worden sind. 
Zum Schluß berichtet der Verfasser über Schwin- 
gungen an Bauwerken, insbesondere an Schorn- 
steinen, Glockentürmen und Wehranlagen, so daß 
ein bautechnisch wertvolles, interessantes Buch 
entsteht, mit dem sich der Verfasser den Dank der 
werktätigen Ingenieure errungen hat. Das Buch 
ist ausführlich geschrieben, der Inhalt durch Er- 
fahrungswerte, Zahlenrechnungen und gute Skiz- 
zen erläutert, kann also zum Studium bestens 
empfohlen werden. 


Dresden. K. Beyer. 523 
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Dr.-Ing. habil. Dr. E. h. WILHELM NUSSELT, 
o. Prof. a. d. Techn. Hochschule München, Tech- 
nische Thermodynamik, I. Grundlagen. 
2, verb. Aufl. (Sammlung Göschen, Bd. 1084.) 143 8. 
m. 71 Abb. Berlin 1943, Verlag Walter de Gruyter & 
Co. Preis geb. 1,62 M. 

Nach ausführlicher Darlegung der beiden Haupt- 
sätze wird die Thermodynamik der Gase und Dämpfe 
etwa in dem Umfang und der Darstellung behandelt, 
die der Einführungsvorlesung in dieses Gebiet an 
Technischen Hochschulen entspricht. Gegenüber der 
1. Auflage wurde der Abschnitt über den II. Haupt- 
satz erheblich gekürzt und noch anschaulicher dar- 
gestellt. Neu angefügt wurde ein Abschnitt über 
die strömende Bewegung der Gase und Dämpfe. 

Die Darstellung zeichnet sich durch eine ebenso 
knappe wie ungemein klare Sprache aus und bietet 
dem Studenten neben der Vorlesung eine hervor- 
ragende Zusammenfassung dieses Wissensgebietes. 
Für Kriegsteilnehmer wird das Bändchen ein be- 
sonders willkommenes Geschenk sein. 

Dresden. Pauer. 506 


MAX LANDOLT, Prof. u. Direktor des Techni- 
kums des Kantons Zürich in Winterthur, Größe. 
Maßzahl und Einheit. 8 S. Zürich 1943, 
Rascher Verlag. Preis fres. 5,8. 

Das Rechnen mit physikalischen Größen, d. h. 
mit dimensionsbehafteten Zahlen, ist für den Phy- 
siker und Techniker ständig wiederkehrende Not- 
wendigkeit. Verf. will in seinem Buche 'die auf 
diesem Gebiet auftretenden Probleme aufzeigen 
und lösen. Er gibt zunächst im ersten Teil seiner 
Schrift eine klare und gründliche Anleitung zum 
Operieren mit diesen Größen. An zahlreichen ein- 
fachen Beispielen zeigt er, welche Rechengesetze 
man auf die physikalischen Größen anzuwenden 
hat, er stellt die Gleichungen auf, denen sie ge- 
nügen, bringt die Umrechnungsmethoden in andere 
Maßeinheiten usw. Zum Schluß folgen dann noch 


eine Reihe vollständig durchgerechneter Übungs- 


beispiele. Der zweite, umfangreichere und bedeut- 
samere Teil bringt die mathematische Begründung 
jenes Größenkalküls. Verf. zeigt mit Zuhilfenahme 
der Gruppentheorie und unter wiederholtem Hin- 
weis auf die physikalische Erfahrung, daß die Auf- 
fassung der physikalischen Größen als wirklicher 
Produkte aus Zahlenwert und Maßeinheit ihre Be- 
rechtigung hat, und daß man zumindest im Be- 
reich der ersten sechs Grundoperationen die für die 
rationalen Zahlen geltenden Gesetze auch auf die 
physikalischen Größen anwenden darf. — Die 
beiden Teile der Schrift sind abgeschlossene Ganze 
und können unabhängig voneinander gelesen 
werden. 
Berlin. Mosch. 516 
Dr. WOLFGANG FINKELNBURG, Prof. a. d. 
Reichsuniversität Straßburg, Physik (Wissen- 
schaftskunde in Einzeldarstellungen, herausgegeb. 
von Dr. Wolfgang Strauß u. Hans Ferdinand 
Schulz, Heft 3). 47 S. Freiburg i. Br. 1942, Hans 
Ferdinand Schulz Verlag. Preis kart. 1,80 M. 
Verf. will in dieser Schrift dem der Physik und 
ihren großen Erfolgen Interesse entgegenbringenden 
Nichtfachmann einen Überblick über Art und 
Inhalt dieser Wissenschaft geben. Er schildert zu- 
nächst Wesen und Bedeutung der Physik im 
Rahmen der Naturwissenschaften, erläutert dann 
die Einteilung in experimentelle und theoretische 
Physik und legt die Methoden der physikalischen 
Forschung dar. Sodann geht er zur Behandlung 
der Einzelgebiete über. Nach Besprechung der 
wichtigsten Begriffe und Sätze aus Mechanik. 
Akustik, Optik und Elektrizitätslehre, aus dem Ge- 
biete der „Makrophysik“ also, versucht er, dem 


l.aien auch eine Vorstellung von dem jüngsten Ent- 
wicklungsstadium der Physik zu verschaffen, der 
„Mikrophysik“, die ja noch in keiner Weise ein 
fertiges, abgeschlossenes Gebilde darsteilt und 
nicht nur dem Laien, sondern auch dem Forscher 
manche stark befremdende Schwierigkeit bietet. 
Dieses Kapitel, das sich u. a. mit dem Kor- 
puskularcharakter des Lichtes und der Wellennatur 
der Materie, mit der Heisenbergschen Ungenauig- 
keitsrelation und der Quantenmechanik beschäftigt. 
kann natürlich nur eine Ahnung von der Größe und 
Wichtigkeit der physikalischen Forschung unserer 
Tage geben, damit aber auch einen Anreiz, sich 
mit diesen Dingen näher zu beschäftigen. Nach 
einer knappen Schilderung der Grenzgebiete der 
Physik (physikalische Chemie, Astrophysik, Geo- 
physik, Biophysik, medizinische Physik), widmet 
er dem aus dem heutigen Stande der Forschung 
sich ergebenden naturwissenschaftlichen Weltbild 
eine kurze Darstellung. Die Schlußkapitel bringen 
Ausführungen über den Beruf des Physikers und 
über das physikalische Schrifttum. Wie man aus 
diesen kurzen Angaben ersieht, ist das anschau- 
lich und flüssig geschriebene Buch von großer Viel- 
seitigkeit. Insbesondere für den angehenden Phy- 
siker, der von der Schule her schon einige Kennt: 
nisse über Physik mitbringt, und der sich einen 
Begriff von Umfang und Möglichkeiten des ge- 
wählten Berufs machen möchte, dürfte das Buch 
einen recht guten Leitfaden abgeben. 


Berlin. Mosch. 517 


Dr. WILHELM MAGNUS, Prof. a. d. Techn. 
Hochschule Berlin, u. Dr. FRITZ OBERHETTINGER, 
Formeln und Sätze für die speziellen 
Funktionen der mathematischenPhy- 
sik. (Die Grundlehren der mathematischen Wissen- 
schaft in Einzeldarstellungen, Bd. LI.) VII + 
173 S. Berlin 1943, Springer-Verlag. Preis geb. 
13,20 M. 

„Die zunehmende Verwendung mathematischer 
Hilfsmittel in der physikalischen und technischen 
Literatur macht ein ständiges Nachschlagen in um- 
fangreichen, mitunter schwer zugänglichen Werken 
der mathematischen Fachliteratur und in zahl- 
reichen Einzelarbeiten notwendig.“ Mancher wirıl 
sich daher schon für seinen eigenen Gebrauch ein 
Verzeichnis der immer wiederkehrenden Formeln 
angelegt haben. Diese Arbeit erspart das vor- 
liegende Buch. Es gibt eine sorgfältig ausgewählte, 
formelmäßige Zusammenstellung der Eigenschaften 
einer Reihe von speziellen Funktionen, die bei der 
geschickten Anordnung und bei dem übersicht- 
lichen Druck ein leichtes Auffinden etwa gesuchter 
Formeln ermöglicht. 

Man findet in den ersten Abschnitten Sätze und 
Formeln für die Gammafunktionen, die hypergeome- 
trischen Funktionen und für die Zylinderfunktionen; 
hier wird auch kurz das Wichtigste über Mathieusche 
Funktionen gegeben. In einem weiteren Kapitel 
sind die Formeln für die Kugelfunktionen nebst 
Kegel- und Thorusfunktionen sowie die Funktionen 
von Gegenbauer zusammengestellt. Als Beispiel 
für weitere orthogonale Polynome findet man die 
Formeln für die Polynome von Tschebyscheff, Her- 
mite, Jacobi und Laguerre. Das sechste Kapite! 
bringt dann die konfluente hypergeometrische Funk- 
tion und ihre Spezialfälle, wie die Laguerreschen 
Funktionen, das Fehlerintegral, die Fresnelschen 
Integrale, Integrallogarithmus, Exponentialintegral, 
Integralsinus und -kosinus usw. Es folgen die 
Thetafunktionen, elliptische Funktionen und Inte- 
erale, weiter die Integraltransformationen und Inte- 
eralumkehrungen mit Tabellen, vor allem für die 
Fourier- und Laplacetransformation, ferner ver- 
schiedene Beispiele für Integralgleichungen erster 
Art. Im letzten Kapitel werden Differential- 
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operatoren für verschieden orthogonale Koordi- 
naten aufgeführt, Beispiele für die Trennung der Ver- 
änderlichen gegeben und kurz auf lineare Differen- 
tialgleichungen zweiter Ordnung eingegangen. Eine 
Verzeichnung der Abkürzungen und Funktions- 
symbole erleichtert die Benutzung, ein Literatur- 
verzeichnis ermöglicht, wenn nötig, ein leichtes 
Aufsuchen der Ableitungen und Beweise. Tabellen 
für einen großen Teil der hier behandelten Funk- 
tionen finden sich in den Tafeln von Jahnke-Emde, 
die ja auch einige der wichtigsten Formeln bringen. 
Das vorliegende Buch gibt eine weit ausführlichere 
Zusammenstellung dieser Formeln. Beide Werke 
ergänzen sich so in sehr nützlicher Weise. 

Die Verf. haben mit diesem Buch ein Werk ge- 
schaffen, das wirklich fehlte. Der große Kreis von 
Mathematikern, Physikern und Technikern, denen 
bei ihrer Arbeit die hier behandelten Funktionen 
vorkommen, wird ihnen für die mühevolle Arbeit, 
die die geschickte Zusammenstellung und die Prü- 
fung des umfangreichen Formelmateriales erfordert 
hat, sicher dankbar sein. 

Dresden. Willers. 512 

Profi. Dr. RUDOLF ROTHE, Höhere Mathe- 
mathik für Mathematiker, Physiker 
und Ingenieure, Teil IV: Übungsauf- 
gaben mit Lösungen, Formelsamm- 
lung. 7. Heft, Zusammenstellung der 
Formeln und Lehrsätze, unter Mitwirkung 
von Fachstudienrat Oskar Degosang u. Dr.-Ing. 
Gerhard Dobrack. (Teubners Mathematische 
Leitfäden, Bd. 43.) VI+ 120 S. m. 74 Abb. Leipzig 
und Berlin 1943, Verlag B. G. Teubner. Preis 
kart. 4 M. 

Mit dem vorliegenden 7. Heft des vierten Teil- 
bandes ist das weitverbreitete Werk von R. Rothe, 
dessen erster Band in erster Auflage 1927 erschien, 
zum Abschluß gekommen. Der Verf. hat die Voll- 
endung seines Werkes leider nicht mehr erlebt, doch 
hat ihm noch das druckfertige Manuskript dieses 
letzten Heftes vorgelegen. Es bringt in systema- 
tischer Zusammenstellung die wichtigsten Defini- 
tionen, Sätze und Formeln unter Hinweis auf die 
entsprechenden Paragraphen des Buches, kann aber 
auch unabhängig davon sehr gut als Formelsamm- 
lung gebraucht werden. 

Das ausgezeichnete Werk, dessen erste Bände 
schon eine Reihe von Auflagen erlebt haben, 
wird auch weiterhin insbesondere von Studieren- 
den neben den Vorlesungen gern gebraucht werden. 

Dresden. Willers. 505 


Dr. J. HORN, em. o. Prof. a. d. Techn. Hoch- 
schule Darmstadt, Gewöhnliche Differen- 
tialgleichungen, 4., fast unveränderte Auf- 
lage. (Göschens Lehrbücherei, 1. Gruppe, Reine 
und angewandte Mathematik, Bd. 10.) VIII + 1968. 
m. 4 Fig. Berlin 1943, Verlag Walter de Gruyter & 
Co. Preis geb. 10 M. 

Wie die 1937 erschienene dritte Auflage ist auch 
die jetzt vorliegende vierte ein fast unveränderter 
anastatischer Neudruck der 1927 erschienenen 
zweiten. Es erübrigt sich, der Neuauflage dieses 
viel gebrauchten einführenden Lehrbuches emp- 
fehlende Worte mit auf den Weg zu geben. 


Dresden. Willers. 513 


Dr. A. TIMPE, o. Prof. a. d. Techn. Hochschule 
Berlin, Höhere Mathematik 1, II. Auto- 
graphierte Vorlesungen. 84 u. 130 S. Stuttgart und 
Berlin 1942, Verlag W. Kohlhammer. Preis brosch. 
Teil I 3 M., Teil II 3,50 M. 

Inhalt: Höhere Mathematik I. 1. Funktionenlehre. 
2. Grenzwert und Stetigkeit. 3. Differentialrechnung. 
4. Mittelwertsatz, unbestimmte Formen, Taylorsche 


Formel, Kurvendiskussion. 5. Integralrechnung. 
Teil U. 1. Funktionen von mehreren Veränder- 
lichen. 2. Unentwickelte Funktionen, Differential- 
geometrie der ebenen Kurven. 3. Linien- und Ge- 
bietsintegrale. 4. Unendliche Reihen. 5. Komplexe 
Veränderliche und Funktionen. 

Autographierte Vorlesungen können immer nur 
ein Notbehelf sein, um den empfindlichsten Mangel 
an Lehrbüchern zu mildern. Darnach darf an eine 
solche Vorlesung natürlich nicht derselbe strenge 
Maßstab gelegt werden wie an ein Lehrbuch. 
Trotzdem m man in dem vorliegenden Fall 
ein Auge zudrücken, zum mindesten, was 
das Äußere der Timpegchen Vorlesungen betrifft: 
Auf 214 Textseiten finden sich rund 50 (!) hand- 
schriftliche Verbesserungen, der Text ist an ent- 
scheidenden Stellen (Grenzen der Doppelintegrale 
usf.) schlecht lesbar, die Zeichnungen sind nicht 
immer gut geraten. Auch inhaltlich kann man 
nicht mit allem einverstanden sein, z. B. wenn der 
Satz, daß eine algebraische Gleichung n-ten Grades 
höchstens nWurzeln hat, als Fundamentalsatz der 
Algebra bezeichnet wird oder wenn behauptet 
wird, daß die trigonometrische Reihe Fb, cos nz 
mit positivem, monoton zu null abnehmendem b, 
für alle z konvergiere. 

Den Hörern des Verfassers mag das Heft an 
Stelle eines eigenen Manuskriptes und als Ergän- 
zung des Gehörten vielleicht gute Dienste leisten, 
einem Außenstehenden kann es jedoch kaum emp- 
fohlen werden. 

Freiberg. Grüß. 504 

Prof. Dr. FR. A. WILLERS, Mathematische 
Instrumente. 270 S.m. 199 Bildern. München 
u. 1943, Verlag R. Oldenbourg. Preis geb. 
14M. 

Das Buch entstand durch Überarbeitung und Er- 
änzung verschiedener Aufsätze des Verfassers, 
die im Laufe der letzten Jahre im Archiv für 
Technisches Messen erschienen sind. Der Ver- 
lag, von dem der Wunsch hierzu ausgegangen ist, 
hat sich durch diese zusammenhängende und er- 
weiterte Herausgabe der Aufsätze in einem sehr 
gut ausgestatteten Buch ein großes Verdienst er- 
worben. 

Das Buch ist in sechs Abschnitte eingeteilt, in 
denen Rechenschieber, Rechenmaschinen, Koordi- 
natographen sowie Kurvenmesser und Differen- 
tiatoren, Planimeter, harmonische Analysatoren 
und Integraphen sowie Integratoren behandelt 
werden. Ein umfangreiches Verzeichnis des in 
Frage kommenden Schrifttums sowie ein Namen- 
und Sachverzeichnis sind beigefügt. 

Der Abschnitt über Rechenschieber umfaßt nur 
zwei Seiten. Mit Rücksicht auf die große prak- 
tische Bedeutung des einfachen Rechenschiebers 
und der weiteren rechenschieberartigen Rechen- 
vorrichtungen wäre eine Erweiterung dieses Ab- 
schnittes in einer späteren Auflage zu begrüßen. 

Der Abschnitt über Rechenmaschinen ist ein- 
geteilt in Addiermaschinen, erweiterte Additions- 
maschinen und eigentliche Multiplikationsmaschi- 
nen. Im Zusammenhang: damit werden einige Re- 
geln für das Maschinenrechnen mitgeteilt. Bei den 
reinen Addiermaschinen wird zwischen einstufig 
arbeitenden mit Volltastatur, zweistufig arbeiten- 
den mit Zehnertastatur und nach dem Speicher- 
verfahren arbeitenden Maschinen unterschieden. 
Die zum Multiplizieren und Dividieren verwend- 
baren erweiterten Additionsmaschinen sind ein- 
geteilt in Sprossenrad-, Staffelwalzen-, Propor- 
tionalhebel- und Schaltklinkenmaschinen. Bei allen 
Bauarten werden außer ihrem Aufbau die ihnen 
eigenen Einrichtungen und ihre Verwendungsmög- 
lichkeit zur selbsttätigen Multiplikation und Di- 
vision beschrieben. 
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Im dritten Abschnitt werden Koordinatographen 
zum Auftragen nach rechtwinkligen Koordinaten 
und nach Polarkoordinaten, Kurvenmesser und 
Differentiatoren behandelt; bei den letzteren han- 
delt es sich in der Hauptsache um Vorrichtungen 
zum Zeichnen der Tangente in einem Punkt an 
eine Kurve. 

Der Abschnitt über Planimeter beginnt mit 
einem kurzen geschichtlichen Überblick und be- 
faßt sich mit den Vorrichtungen für die Flächen- 
bestimmung durch Zerlegen und Verwandeln, mit 
der Theorie der Umfahrungsplanimeter, den Inte- 
griervorrichtungen den Polarplanimetern, den ver- 
schiedenen Linearplanimetern, den Linear-Potenz- 
planimetern, den Radialplanimetern, dem Schneide- 
planimeter und den Integrimetern. Anhangsweise 
werden die Fehlerursachen bei Messungen mit 
einer Integrierrolle behandelt. 

Der den harmonischen Analysatoren gewidmete 
Abschnitt ist eingeteilt in Bestimmung der Fou- 
rier-Koeffizienten mit dem Planimeter, harmo- 
nische Analysatoren, harmonische Analyse mittels 
gleichabständiger Ördinaten und das Stieltjes- 
Planimeter. 

In dem letzten Abschnitt werden Grundintegra- 
phen und Differentiatoren, Anwendungen der Inte- 
graphen, allgemeine Integraphen sowie Integrier- 
maschinen behandelt. 

Wie diese kurze Übersicht zeigt, ist der Inhalt 
des Buches überaus reichhaltig, so daß jeder, der 
sich für mathematische Instrumente interessiert, 
das ihn Betreffende darin finden wird. Gute Ab- 
bildungen und klare, schematisch gehaltene Dar- 
stellungen sind reichlich vorhanden. 

Da die Zahl der Bücher über mathematische In- 
strumente nicht groß ist, und das vorliegende 
Buch nach Inhalt und Behandlung des Stoffes 
bestens empfohlen werden kann, so wird es sicher 
viele geben, die das Buch mit Nutzen zur Hand 
nehmen. 

Stuttgart. P. Werkmeister. 527 

Prof. Dr. HERMANN SCHUBERT, Mathe- 
matische Mußestunden, eine Samm- 
lung von Geduldspielen, Kunst- 
stücken und Unterhaltungsaufgaben 
mathematischer Natur. Neubearb. von 
Prof. Dr. Fitting. 10. vermehrte Auflage. 271 S. 
Berlin 1943, Verlag Walter de Gruyter & Co. Preis 
geb. 4,50 M. 

Das beliebte Schubertsche Buch, das in der 
Neubearbeitung eines so ausgezeichneten Kenners 
wie F. Fitting außerordentlich gewonnen 
hatte, erscheint hier nochmals erweitert und be- 
reichert in zehnter Auflage, 271 Seiten stark. Wir 
verdanken F. Fitting interessante und wichtige 
Forschungen auf dem Gebiet der magischen Qua- 
drate und können es nur lebhaft begrüßen, daß er 
einige seiner neueren Ergebnisse in diese neue 
Auflage hineingearbeitet hat. Auch das Kapitel 
über Rösselsprünge bietet in der vorliegenden Auf- 
lage manches Neue, so z. B. ein einfaches Ver- 
fahren, um auf dem Schachbrett mit 36 Feldern 
alle möglichen Rösselsprünge zu finden. Das Buch 
wird besonders unsern Frontsoldaten und Ver- 
wundeten eine sehr willkommene Lektüre sein. 
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Dr. WL. K. HRISTOW, Chef d. Astron. Dienstes 
d. Militärgeogr. Instituts zu Sofia, Die Gauß- 
Krügerschen Koordinaten auf dem 
Ellipsoid. VI+80S. m. 9 Abb. Leipzig und 
Berlin 1943, Verlag B. G. Teubner. Preis kart. 
5M. 


Die vonK.F.Gauß bei seiner Landesvermessung. 


von Hannover benutzten, vom Ellipsoid ausgehen- 
den ebenen Koordinaten, durch die die winkel- 


treue zylindrische Abbildung in transversaler Lage 
bestimmt ist, wurden insbesondere vonL. Krüger 
einer eingehenden: Bearbeitung unterzogen; sie 
werden deshalb seither als Gauß-Krügersche Ko- 
ordinaten bezeichnet. Nachdem das betreffende 
Werk von L. Krüger!) und auch ein früher er- 
schienenes, denselben Stoff behandelndes von 
O0. Schreiber?) schon längere Zeit vergriffen 
sind, hat der Verfasser mit Benutzung eigener, 
bereits veröffentlichter Vorarbeiten in dankens- 
werter Weise die Aufgabe übernommen, eine um- 
fassende und gründliche Darstellung dieser Koor- 
dinaten zu geben. 

Der Stoff ist eingeteilt in einen theoretischen Teil 
und in einen praktischen Teil. Der letztere Teil 
enthält außer einer Zusammenstellung aller For- 
meln eine Reihe von Tafeln, denen das Besselsche 
Erdellipsoid zugrunde liegt. Am Schluß werden 
für alle in Frage kommenden Aufgaben durch- 
gerechnete Zahlenbeispiele mitgeteilt. 

Mit Rücksicht "auf die nach Einführung der Me- 
ridianstreifen große Bedeutung der Gauß-Krüger- 
schen Koordinaten und der durch sie bestimmten 
Abbildung hat sich der Verfasser durch seine Ar- 
beit ein großes Verdienst erworben. 
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1) Konforme Abbilduug des Erdellipsoids in der Ebene. 
Potsdam 1912. « ps in der Ebene 


2) Theorie der Projektionsmethode der Hannoverschen 
Landesvermessung. Hannover 1866. 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


F. KOHLRAUSCH, Praktische Physik 
zum Gebrauch für Unterricht, For- 
schung und Technik, Bd. 1 u. Bd. 2. Her- 
ausgegeben von F. Henning, Direktor bei der 
Physikalisch-Technischen Reichsanstalt. 18., neu- 
bearb. Aufl. VII+535 S. m. 298 Abb. u. VI+ 
578 S. m. 328 Abb. Leipzig u. Berlin 1943, Verlag 
B. G. Teubner. Preis zus. 35 M. 


Dr. phil. habil. U. STILLE, Dozent a. d. Techn. 
Hochschule Braunschweig, Umrechnungs- 
tafel für die Zahlenwerte elektri- 
scherund magnetischer Größen. Braun- 
m 1943, Verlag Friedr. Vieweg & Sohn. Preis 
250 


Dr. phil. habil. H. GEBELEIN, Baurat a. d. Staatl. 
Ing.-Schule Eßlingen u. Dozent an d. Techn. Hoch- 
schule Stuttgart, Zahl und Wirklichkeit, 
Grundzüge einer mathematischen 
Statistik. Mit einem Geleitwort von Geh. Reg.- 
Rat Prof. Dr. E Wagemann, Präsident des 
Deutschen Instituts für Wirtschaftsforschung. 
+425 S. Leipzig 1943, Verlag Quelle & Meyer. 
Preis geb. 11 M. 


Dr. phil. GEORG PRANGE, weiland Prof. a. d. 
Techn. Hochschule Hannover, Vorlesungen 
über Integral- und Differential- 
rechnung, I. Bd. Funktionen einer re- 
ellen Veränderlichen. Herausgegeben von 
Dr. phil. Werner v. Koppenfels, Prof. a. 
d. Deutschen Technischen Hochschule Brünn. XI 
+ 486 S. m. 140 Abb. Berlin 1943, Springer-Verlag. 
Preis geb. 22,50 M. 


Prof. Dr. BERNHARD BAULE, Die Mathe- 
matik des Naturforschers und In- 
senieurs, Bd. II, Analytische Geo- 
metrie. 80 8. m. 89 Abb. Leipzig 1943, Verlag 
S. Hirzel. Preis brosch. 4,40 M. Del. Bd. IV, 
Gewöhnliche Differentialgleichun- 
ren. 102 S. m. 41 Abb. Leipzig 1943, Verlag 
S. Hirzel. Preis brosch. 5,80 M. 
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ZUSCHRIFTEN AN DEN HERAUSGEBER 


F. Reutter: Eine Anwendung des absoluten 
Parallelismus auf die Schalentheorie. Z. angew. Math. 
Mech. Bd. 22 (1942), S. 87 bis 98. 

In obigem Aufsatz werden am Schluß die Fun- 
damentalgrößen 2. Ordnung L, M, N für eine de- 
formierte Schalenmittelfläche abgeleitet. Für M 
wird die Formel angegeben 


Hierin bedeuten v, die Verschiebungskomponen- 
ten in Richtung von 


Zu Iv_ 
VE 
die partiellen Ableitungen nach den Parametern 
u und v. Es soll nun gezeigt werden, wie Gl. (1) 
in eine für weitere Rechnungen sehr übersichtliche 
Form gebracht werden kann. 
Zunächst formen wir in .(1) die symbolischen 
Operatoren um, indem wir schreiben: 


‚v und Yu, Un 


Nach entsprechender Zusammenfassung kommt aus (3) 

—: —u+Eu-— | —— 


L 
2 — 
VE 


(4). 


vs 


Entsprechend der Voraussetzung, daß die Para- 
meterkurven der undeformierten Schale Krümmungs- 
liniien (’=M=0) sein sollen, erhalten’ wir aus 
den Mainardischen Gleichungen !) 
a... 
die weitere Form 

N 


dulgyG) e2eyG\E 
Aus (4) wird dann in Verbindung mit 
LE D) u 
VE ' EYE dv\yE)’ 
(77). 
vs " du\yG 


Schließlich erhält man die sehr bequeme Darstellung 
von M in der übersichtlichen Form: 


Man erkennt hierin sehr schön den Zusammenhang 
von M mit den alten Fundamentalgrößen und den 
Verschiebungskomponenten. 


Karlsruhe, z. Zt.imFelde. J. FadleVDI. 511 


1) Man vgl: z. B: Kommerell, Theorie der Raum- 
ze. .- krummen Flächen. Berlin und Leipzig 1931, 


U. T. Bödewadt: Die Drehströmung über 
festem Grunde. Z. angew. Math. Mech. Bd. %0 
(1940) S. 241—253. 

Die in dieser Arbeit auf S. 242 oben herange- 
zogene Vernachlässigung der Grenzschichttheorie, 
welche das Ergebnis der ganzen Arbeit auf hin- 
reichend große Reynoldssche Zahlen be- 
schränkt, ist vermeidbar. Macht man für den 
Druck p den Ansatz 


p=0 


mit der noch zu bestimmenden Funktion P (£), so 
kommt man genau so zum Ziele, und die weiteren 
Ergebnisse von Bödewadt bleiben unverändert 
bestehen. Die Bödewadtsche Lösung ist also 
eine strenge Lösung der vollständigen Navier- 
Stokesschen Gleichungen und damit für alle 
Reynoldsschen Zahlen gültig. 

Bei dieser Gelegenheit sei auf folgenden Druck- 
fehler hingewiesen. Die in den Gln. (1.2.1) und 
(1.2.2) der obengenannten Arbeit auf den rechten 
Seiten angegebenen Reibungaglieder sind unvoll- 
ständig. Gerechnet wurde jedoch mit den rich- 
tigen Gleichungen. 

Göttingen. H. Görtler. 

Die Zuschrift hat Herrn U. T. Bödewadt vorge- 
legen. 535 
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